FA NOVA R LTE 


NOI raven.) 


psa kanenh eh 
PATE CURD 
eM 


ats Aa et 


balla 


errr rs 


Cea 


pnts ate 


gail 


Is \eaip e Wie ied wal 


ig AST 
wonrhw 


ye Ars 
9 bauriet 


re 
ea A Oat 
Ava nathect ad 


Hereadrnavay) 
ails 


AA nings hey 


a 
fei cigs) 
yea Hf Y, 


He iReas 
Rolteh \elater 
a eal 


o « Area 
a helen 


elf») 


eM by 
hi i 
ny 


ite i 
Bsa 
ieaaint ite 

Mia 


pee 


Gah esifous ea 
Si atenieating nhs 


eg 
& Nfeke ras (ened 
oes Fehr ne 


‘i 
ee ait! 
ipsa? ah 


Ap ahek el 


feotec@ atest 


eh sini 


aay Uaneiaus 
ca 


i 
Mh peat ae 
ye sakes 


ah ons gee 
Tes rena ee Hge te 
Maiieceaa (aie egeipempe ae 2 
S pawilanone rstraeengyay? 
pate hede,sa rey 
ematec rely oh 4 
echo adie 


sinmer aah) 
Apcicapin me ey Alp 
hpeqanavansiione hamapneat 7) 
qed sh ea st eps neage aie = 
Hie piece 


nora ha 
ers heey Aah afeisanene at 
Qatar has yee aaa pater? 
ig FOE AO AN: Hi ype ys 
iastniae ia i ie aay gay ae ver 
ell, 


ciatinekibeined 


bob eke Raina 
Ayer aire! 


VAteine 


ped ieg ven ie 
ee aes 


ae hat “ 
S Kehaker) 
Phe if 
Blas ee 
NF is ele Ba 
eM a4 | lane 

(ep pana 

jieisechige | 


as sicietel 
iraiep haste 
agree 


Beis Pde 

vijuinse to shale! fifa) 
igh peal caer ai! 
abe d dehy 

i paren 


gs 


‘ | 
Creu 
nay 


i ee 


i 2a Nua halenthen 
noel th eae eae uacads 
erga ices aka 
a 


perl ies 
arf cbs iuihoitele 


ui ra ren 
ieee 4a 
Ais 

a 


Rote 
Waibaikaiy Bukit 


intRy VARIN ee 
DIAweNiatey 
BIEL oy 


pyeeTyast date gent alle} 


vaheaitia: wa bien 


ate et" 
pants th 
#4 


A 
Gasuitaid oi) 
LRN VAbaw eS p ED Vase 


Fulignee wey the 
ovale eajent 


praia 


Twine sk 


Sain 
shapes a ak 


Avie 


lieaarkty © 
HLT NOREEN SR OESEE 


Hannti er 


iueses 
pe 


Weert 


ay Ge BR § 


ned Hm 


pata node oR AME 


fife sala SI Wises 
Rea Gove seameneaciciies era h 


ene wives 


Paneer on YL wR ere Ms 


i eal sboceebucsyane 


sem Ri 


fy wabeitedalbs he) 
rf a 


rita hee ea 
ie dst! 


aauelh {ie a wih rsh Ps 
cee aawaett sna Biya eee me 4 ti ie 


Pia re vibts 
peta ies len ae? oy Ti Biwi Gowers asteny 


afc ee “yailautie 


fabdaree hia et ae 
i iomaagit nits 
ene pe 
eave bartell inane! Lit are ning 
erate a 
shisues 
piacere 


yelheveibee 

grdeurh apcely a 

Vea tare * 
sh we 


ate est 
4 i 


Cat 
tie Bee 


Fe ai ee 


ae 
wrmiy sed ota 
seebeedeie 
id he bate 
Trl gies 
Lael 
fababestenedetahiel 


ads Hees 


verk ait alt ae te 
pee ee <i Seancaaee ores 
i 


ite it fiona a) 


agelh Mi 
hah. alhiacadeck ee 
Ziaipsn cnet 


eee hae 
eal 
a i nee Hi 
< eat 
ia eile 


page Vabae 
Kina ie velivioelt oe 
‘there ay he male bailey ae lack (esisane 
Fevogohen, pbetruea hein 
é Tae nse satin c 
leds geod Ratha 
hig one hs 
t) nf 
Me ci reetetteniaic et 
tbls pale fates 
i ; “hf tide scesyeicrteetaks 
i eehyn ste ened fi 


ek 
an) ibcheds 


he! Shae 
wise nip sepa 


Vnlhshd abe hee 
ch 


ee aa) ahd ch sts 
sb hall aes ot 
lish ghar oh th 


a 

Aner sid 
nph ad dh ee 8? 

as} fae AY aie dhe Ry et ete 

ait W i i 

a ahaa tidak iphins} Betas 


Ny 
i Whee 
Sisiayueeeeaete aVaire oy 
ebrirs aye een Ad 
fl)xd wie leo f: ie 
oH Maina Tero Hd a 
Bokeh aaivah salle ecient 
iS re 


5-61 ll sk ois ainsi n Se rBAaedr eb 1 4% 
igrraneieel ces conte 
4 


Oke fe 
See 


Ai wah os 
igatt ee 


Nabkapste Vardar eibattaN 


ine rat =e BS 
iv OP 


seep ey 
in pai teeeng§ie pine 
pethicthe eae lately 
caf ey Snibe a ar PPS W 
Higthrnve M2) mh 
couenae 


vi nih SS 
ve Baie 


ren 
sient mor pinne 
Fig 08 feats 


ARCA ELE. 9, 

ae Need et nee a h 

ahh alt ri 
mh rae Meare aibatiy hs tal eon Aree 
friend easitene he! eke ER Sibu abate 

Nahe y és (pA Warne fre he! 
7 ah iSite Reps ees ¥ SP 
wehateb hace 


Zr Webi 
dbanan dalled 2 

Pit aaa 
mee vt abana 


asia 
Ne Saidct 
ee 


iad Larne 
eon Eeibas aly! ia 
ie Lit gh drem ck alte bh t fe. ured 
basetansmsmetenaari eh 


ike ing wails ht 


apeavesy 
Ry 


eapareaae i 

idata teks parker weieee 

pie kavcne mat a 

aire ieanhh juceeceneryy 

foeansit™ noe? 
serene escent 

Terabe seletou weakel 

awl eaeprere ys ye 

makina aoe 

sf rnae tenses 

pyawnad sink 


Sgathetet 2 
ee pene be « 
es 


Reka 
ewe sea ORE 


Benes 
ieee he 8x 


Bau SA ing elt eee st 
poi Warceateah 9! 
sah ww 


@awicbh ieee oe 
is bese 
caer neeecity 
ra he 


Seite aiea 
felons 
Sateen 


mS a 
nich Sern oy 


Ae 

toys au a> b all gnats 

es sigan desea Fert 
mene 


Par nas 


Pere eer 
i yanea dig wets 


Tie ahve eo egies 


auncrsa 
Act's iment ah sy 


Ex wenis 
ANINTRSUTATIS 
AUSERTAEDSIS 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Alberta Library 


https://archive.org/details/Fung1983_0 


THE UNIVERSITY OF ALBERTA 


AN EXPLORATORY STUDY OF STORMWATER DROPSHAFTS 


Pi by 
f i. \ 
{ 4 


fi 
Lu) LARRY SIU-KUEN FUNG 


AATHES LS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 


OF MASTER OF “SCIENCE 


THE DEPARTMENT OF CIVIL ENGINEERING 


EDMONTON, ALBERTA 


FALL 1983 


Abstract 

This thesis is a result of part 1 of a 2 part research 
project initiated by the City of Edmonton Water and 
Sanitation Department. The objective of the first part of 
the study is to provide a theoretical evaluation of the two 
types of drop structures used in the City, and to provide 
recommendations for the improvement of their design 
practice, as well as further studies (part 2) to answer 
outstanding questions. 

These structures are called the shaft drop structures 
and the drill drop structures primarily due the different 
methods of construction. They are analysed by dividing the 
structures into identifiable hydraulic components. 
Discussions and analyses include the hydraulics of various 
elements of the structures, the possible flow types and 
controls, balance sizing of the components, the aeration 
problem, and the surcharge conditions. Wherever possible, 
theoretical derivations are verified with published data or 
with some new experiments in the laboratory. 

Based on the analysis, criteria for evaluating the 
hydraulic performance of these structures are provided. 
Modifications and improvements are also suggested. Further 
experimental studies are recommended to tie-in with the 
theory and to develop adequate design procedures for these 


Structures. 
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1. Introduction 


1.1 Background 

In a municipal sewer system, drop structures are used 
to disharge wastewater from a higher elevation to a lower 
elevation. They may be required eBy a variety of reasons. 
Typically, for development in a high relief area, drops are 
required to absorbed extra head losses so as to maintain 
sewers at an acceptable grade. For the discharge of 
stormwater directly into the river located in an’ entrenched 
valley with steep valley walls, drops are required to 
connect the sewer elevation oeohe river elevation. Also, 
ror old municipality, the development of underground storage 
tunnels for the temporary storage of excess stormwater 
during flood peak means the requirement of drops to connect 
the sewer system to the storage tunnels underneath. 

The second example listed above is typical of the 
riverbank area in Edmonton, whereas the last example 
describes the situation in Chicago. Usually, vertical drops 
ae used aque to ithey. limitation of space, and cost of 
construction. In Edmonton, two types of drop structures are 
in use. They are the shaft drop structures and the drill 
drop structures. Their typical geometries are shown in 
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Fig. 1 - Geometry of shaft drop structures 
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Fig. 2 - Geometry of drill drop structures 
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1.2 Outline of Presentation 

The shaft drop structures and the drill drop structures 
are discussed individually. The shaft drop structure can be 
subdivided into three major components: 1) the inlet sewer 
and the shaft entrance; 2) the drop shaft; and 3) the plunge 
pool and the entrance to the outlet sewer. The drill drop 
Structure can also be subdivided into three components: 1) 
the inlet sewer; 2) the entrance chamber; and 3) tthe drill 
drop Pipe. The hydraulics of each component will be 
discussed . 

The various types of flow controls which may exist in 
this structures are presented. Theoretical development of 
the governing equations for these flow types and the 
conditions for their existence is also provided where 
possible. A computer program for the evaluation of the drill 
drop structures is also included bearing in mind that its 
development is based on the assumptions and approximations 
discussed in the text. 

Air entrainment in drop shafts has’ been an active 
research area for the last twenty years. Due to _ the 
complexity of the problem, no comprehensive theoretical 
treatment is available at this time. Furthermore, its 
dependence on viscous and surface tension effects has made 
modelling extremely difficult. Based on recent research 
developments, some modelling criteria are now available. 
However, many problems have yet to be resolved before 


Satisfactory models can be developed. For practical 
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purposes, empirical relations for air demand may be 
developed based on field observations. 

Overloading, in terms of exceeding the design capacity 
and maximum tailwater level, is also discussed. Other 
associated problems, such as erosion due to high velocity, 
Cavitation, impact pressure at the base of drop shaft, and 
the hydrostatic uplift pressure due to high groundwater 
table, are also included. 

Based on the theoretical discussions, some possible 
modifications are suggested for consideration to improve the 
hydraulic efficiency of such structures. Recommendations for 


future experimental studies are also provided. 
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2. Shaft Drop Structures 
2.1 Inlet Sewer 


2.1.1 Discharge Capacity 

The discharge capacity of the inlet sewer may be 
determined with simple culvert hydraulics. In general, the 
discharge and depth in a sewer can be described by the 
following relation, 

yD 4b (0/g)! 7D" HAwSis « Hise) (1) 
where y is the normal depth, Q is the discharge, S$ is the 
Slope of the invert, D is the diameter, and f' is the 
Ericciony factor. pin’ particular, ..the Gcriticaim depth and 
discharge are uniquely related and can be written as, 
OTP mera pears ee (2) 

A culvert may have inlet control, outlet control, or 
Pipe-flow control, depending on the slope, length, entrance 
condition, and the tailwater level. A definition sketch of 
these is shown in Fig. 3. In a normal design situation, the 
discharge capacity and allowable head water-level are 
predetermined. The allowable slope is determined by site 
conditions and other sanitary constraints. If the capacity 
of the vertical shaft is designed so that backup at the 
shaft-entrance will not occur, and if the sewer slope is 
mild, critical depth control at the oulet will always exist. 

One method of solving for the appropriate inlet size is 


to use a trial-and-error procedure. First of all, a trial 
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Supercritical 
flow profiles 


Orifice control, 
inlet submerged 


Weir control, 
inlet unsubmerged 


(A) Inlet control, sewer is short and steep 


Gradually varied flow profile, 
might become uniform flow in the 
upstream reach 


(B) Outlet control when the outlet is unsubmerged; 
sewer slope is mild 


Velocity increase due to 
entrance condition 


(C) Pipe flow control when both the inlet and the 
outlet are submerged 


Figure. 3.- Definition sketch of the inlet flow controls 
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pipe size is selected and the critical depth for the design 
discharge is computed. Then, the flow depths at various 
locations of the pipe are determined by a gradually varied 
flow computation, starting at the critical section near the 
outlet to a section just behind the inlet. The headpool 
water-level may be calculated by taking into account the 
entrance head loss. Tf it 1s different from the required 
headpool level, a new pipe size will be chosen and the 
procedure repeated. 

According to Henderson(1966), a culvert will flow full 
when y,/D is greater than 0.9, except near the brink. In 
this case, the discharge capacity may be computed using the 
full-pipe-flow equation with the head equal to the _ head 
difference between the headpool and the center line of the 


sewer outlet. 


2.1.2 Momentum Considerations 
If the inlet flow goes through critical depth near the 
shaft entrance, the discharge in the inlet may be determined 
by the momentum equation between the critical section and 
the end depth section as discussed by Rajaratnam and 
Muralidhar(1964). The equation can be written as: 
Pee eS TO be et Meee Bt (3) 
where P is the total pressure at the critical section, P, is 
the pressure at the end section, W is the weight of fluid 


between the two sections, a is the angle of the pipe invert 


to vthe horizontal’, Py is the boundary friction force, M, is 
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the momentum flux at the end section, and M is the momentum 
flux at the critical section. 

For small sewer slopes and short control volume length, 
the term (Wsina-P,) is small compared to the pressure and 


momentum terms, hence, 


PeS[Ps iseMet- M (4) 
wherein P = yyAcosa, (5) 
Be = Y¥,A,K1C0Sa, (6) 
ySert 
M = *—°B, (a7?) 
g A 
12301 
and M. = -—-B (8) 
g A, 


where y is the depth from the water surface to the centre of 
gravity of the section, A and A, are the cross-sectional 
areas for the respective sections, and B is the momentum 
coefficient. 

By Ssubstmtucmnge (5), (6), (7), and (8) into (4), and 


putting cosa ~ 1 for small invert slope, (4) can be changed 


into: 
O# 1 1 
WeDo hey at ew OD IX ey = Al yD" at b>? a) 
where A= UD? - and Woe £3: y/D_) (10) 
and VeoAys ended = f4(C y/D ) Ca) 


wherein w and 2d are dimensionless coefficients and are 
unique functions of y/D, K; is the pressure coefficient at 
the end section, A is the area, y is the depth of the centre 
of gravity below free surface, y. is the critical depth, ye 


is the end depth, and, c and e stand for the critical 
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section and the end section respectively. For the critical 
discharge, 


Oo? 
ee) 


Pew! 

Vv D 
From the experiments of Rajaratnam and Muralidhar(1964), 
Bey x O).\725+f0r 0,.< y,/D SU 29 and a horizontal pipe. 
ALSO; Kay iS sa, Lunction of So/S, only, and for a pipe slope 
of zero, K, 18S zero. The variations of ve (bn and K, with 
So/S. are shown in Fig. 4 and Fig. 5. If the discharge in 
the inlet sewer is known, the critical depth and the end 
depth at the shaft entrance can be computed. Combining (9) 


and (12) gives: 


Q ane 1 
a all a: a3) 
¥gD<< Ven ive 
where Bes Buy. / Ds 
x Be & 3 
and Be, A,y,/D 


Putting oO: 


Ogi? 2D! *., “wegets: 
peeesK Pp * leiee 
oe [EE 
C1/v, - We) 


Values of Poly Ber Ve , and Ve may be obtained from 


(14) 


Appendix I. 


2.1.3 Shaft Entrance Conditions 
Unless a smooth vertical curve is provided, the lower 
nappe will not follow the inlet boundary. Instead, it will 


leave the invert to strike the vertical shaft wall. A 
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surface roller will form on the upper nappe which entrains 
air'. Part of the entrained air will be carried down the 
shaft giving the downward flow a foamy-white appearance. If 
air is not allowed into the upper air chamber or the rate of 
air inflow is less than that of evacuation, it may be drawn 
in from the upstream manhole. The pressure in the air pocket 
above the upper nappe will then be somewhat subatmospheric. 
If adequate ventilation is designed for and maintained 
during operation, air pressure above the nappe will be 
almost atmospheric and the air supply will not be a problem. 

Without a model analysis, it is difficult to predict as 
to whether the nappe will separate the upper and lower air 
chambers. For the moment, if separation is assumed to be 
possible when the drop shaft size is the same as the inlet 
size, and when the y,/D ratio is greater than 0.5, then 
there can be two main categories of flow: one in which the 
upper and lower air chambers are connected, and the other in 
which they are separated. A definition sketch of the flow at 


the shaft entrance is shown in Fig. 6. 


2.1.3.1 Flow Conditions When the Flow Separates the Upper 
and Lower Air-chambers 

In this situation, the upper ‘air chamber above the 
nappe is vented by the manhole at ground level, while the 
lower chamber is ventilated by the outlet. If the outlet 
remains open and is ventilated, the downward flow will be an 


' Air entrainment by a surface roller due to impingement has 
been summarized by Rao and Kobus(1974), and will be 
discussed in a later section on air entrainment. 
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Figure 6 - Definition sketch of the flow at the shaft 
entrance 
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aerated vertical flow. » Since . the») pressure «in the shaft 
remains approximately atmospheric, the energy balance in the 
shaft is generally assumed to obey the Bernoulli criterion. 
The air entrainment aspects of this type of flow will be 
discussed in a later section. 

If the drop shaft is lengthened such that the expanding 
downward flow eventually fills the cross-section of the 
shaft before it reaches the outlet, due to the air 
entraining capacity of the flow, it is possible that the air 
pocket in the shaft may gradually be evacuated, forming an 
air-water mixture in the shaft. One-dimensional two-phase 
flow in a vertical pipe has been analysed by Wallis(1969). A 
Summary of the derivation of the momentum equation, together 
with a worked example is shown in section 2.2. 

A homogeneous two-phase flow is probably not a steady 
State phenomenon in an_ open dropshaft system. Once 
substantial negative pressure is developed in the shaft, air 
will probably be forced into the shaft through the entrance. 
This is a well known phenomenon which was described by 
Blaisdell(1958,1952) as air-controlled flow. The flow is 
characterized by the formation of a noisy air core. The 
velocity of air will depend on the pressure difference 
between the entrance and the dropshaft, and can differ 
significantly from the fluid velocity. The insufflation of 
air may be intermittent or continuous and cause violent 
vibration of the structure. The method of analysis of this 


flow is uncertain. The two-phase flow theory might apply 
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just below the entrance with the negative pressure relieved 
by the air insufflation which gradually diffuses into the 
flow. The flow may even separate from the pipe wall as it 
accelerates down the shaft, and will thus be ventilated by 
the outlet. Further experimental work is needed to clarify 


these comments. 


2.1.3.2 Flow Conditions When the Flow Does Not Separate the 
Upper and Lower Air-chambers 

This flow type will occur when y,/D is less than 0.5. 
The pressure in the lower air-chamber will be atmospheric 
Since it is ventilated. Because the air-supply in the 
G@rcpshart)isenota problem, the flow in’ the shaft will be a 
free-entraining self-aerated flow. Its analysis is similar 
to the previous case where the lower air chamber is 
ventilated by the outlet and will be discussed in a later 


section. 


2.1.3.3 Possibility of Orifice Flow 

Orifice flow is a quiescent stage of flow with the 
entrance flowing full and approximately level. Inside the 
pipe, the flow goes through a contraction with very little 
to no air entrainment from the inlet, and remains partly 
full throughout the pipe length. The annular air space is 
continuous and open to the outlet. This type of flow is 
highly inefficient because the discharge capacity of the 
pipe is greatly reduced with respect to the available head 


at the entrance. 
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Therefore, the entrance should be designed in such a 
way that orifice control will not govern the flow. Orifice 
flow may be suppressed by one of the following ways: 1) 
enlarging the entrance; 2) providing a smooth inlet; or 3) 
insuring that the length of the vertical pipe is long enough 
so that the expanding flow will be capable of evacuating the 
pare spocKet ~in the. shaft. if sorifice flow does occur. 
Experiments done by Anderson and Dahlin(1975) for the Tunnel 
And Reservoir Project (TARP) in Chicago showed that if the 
entrance to the shaft is a smooth elbow, the nappe will 
always follow the inner pipe boundary, and hence weir flow 
will always occur. 

Although the entrance in the present case is not 
smooth, it iS uncertain whether orifice flow would occur. 
Since the form of flow control is strongly dependent on the 
geometry of the inlet, it is advisable to use a smooth 


transition. 


2.1.3.4 Entrance Hydraulics 
Referring to Fig. 6, if we write the Bernoulli equation 
between section® c and a ‘section in front of the roller, we 
have: 
pies v*/2g + AZ = V¥4/20e% Uy, 75} 
where Ve is the critical velocity, Az is the vertical 
distance between the invert of inlet and the mid-depth of 
the flow section in front of the roller, V is the velocity 


in front of the roller, and H, is the head loss between the 


two sections. 
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The surface roller above the upper nappe will reduce 
the kinetic energy of the flow somewhat. If the velocity of 
the downward flow after impingement is V,, the vertical 
component of the momentum equation for the impingement 


region can be written as: 


pQ:1Vm(1-cosé@) PQ3Vm(1+cosé) (16) 


Q3 ia COS 
Or — = (17) 
Q; t+ cos? 
where Q;,, Q3, and @ are defined in Fig. 6. Also, for the 
impingement region: 

PQV = PQi1Vm + PQ3Vm (18) 
where Vm is the velocity just below the impingement region, 
and V is the velocity of the flow before entering the 
roller. Combining (17) and (18) and simplifying: 

(i cose. 4} 
Ve = a WY (19) 
2 
The momentum equation in the x direction is, 
y = 20) 
hee POV pO. ( 


which may be rewritten as: 


r Q 
Bobet. stan hygstiyes) G21) 
D gD? 


The first term is a dimensionless pressure force term 
which is a unique function of y,/D, and was defined earlier 


asec (see Appendix I). The critical velocity, V,, is equal 


to Q/W,D’. Hence, (21) may be rewritten as: 


gD" Ps Q 
Vsing = —— + 
Q wD? 


Cc 


(22) 


7 Te ere an be anes ru ag eee nadine fee Wt 


iGr 
Abe : cal 
enwber fii egies 6a) atid iors ‘oat 


7 “af 406i mt Sas woe! ia 


2 ) ite, oe y 7 
rapepriqnt ba? wot | wie  wosctemen- ahs 269m 


: ‘ e 
ry Pie, aS m= 1 SV TOG ; 
I p 
2 Th 
P Ja 5 J 
. 
} f } { ier é = 7 yt a0’ 
ie ge 
i i ; a ry Ff ¥ On 
; } iy \ m¥ es Me 
‘ 3 : 13. woLem sgRaG, ¥Is 90 e = Fe 
! 7! i} 4 
i : & «& ai : tHe! ease * 4 
i je i 4 
B A = 
: ke Sne 181) Bie (Tt) polerame 
iY ~ 
oy, 
i ( ABpo aT ) - 
‘ = Soy — 
[ if vi . ——e eg er ed = wa ea 1 
Ms ia! 9 ! ha is 
==3 4 ’ ; 
Vee ; \ j 
+ 


i a 


: ; Mee ent ae | 
Bt iif; 5ert6 - 4 ‘an? it » 95 £ Laue dion f mn m I 
i) 


| reer iy 
- +4 . 
(o¢) aE sv ee ae Le at ; ‘ 
ie i a 
as 
; Nie 
) Te 
. 


ot - eign a aves 

a aes & Nid, tee sig I a ne 

= a EP ae apes izpuhchanale 
ech ae bie . 


earied 26% ons VO) 


er: a? 
7 ee ie hata tingiates . 
ed $ee err, ieee’ 


Noe a 


r a 
- ce 2 * i 
a a af sf 3 


18 


Neglecting the head loss term in (15), Substituting it into 


(22) and rearranging: 


+ 
! Q yD? 
sin@ = (23) 
OF 


2gD*y .? 


gD* Pa" Q 


rs 
pea ez. + Ya 


Now, substitute (23) into (19): 


+ 
1 / Q YD? 
Vere ol ee to 
2 Je Q? 
7: gd. “Azer. “yah ) 
2gD‘y 2 
hk Q 

| si oan pees (24) 

V Cc 2gD‘*w .? 


The value of Az should be obtained from experiments. 
With the lack of model observations, Az may be assumed to 


be: 


1 Di Hh aibehe 
Az = --3| ee (25) 
2 Vx 


where H is the entrance step height, and Lr is the length of 
the roller. Combining (21) and (25), and rearranging: 


1 orD thal taelas 1 
Az = —*ge——_—_______—_ (26) 
2 GDP ee Q ; 
oe 


Q y .D? 


If the discharge is known, the average velocity just below 


the entrance, Vm, can be computed. 
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2.2 Considerations for Flow of Air-water Mixture 


2.2.1 Derivation of Momentum Equation 

One dimensional two-phase flow theory was discussed by 
Wallis(1969). The derivation of the steady state momentum 
equation for an air-water flow in a vertical shaft is as 
follows: 

Considering a section of the fluid in the vertical pipe 
as shown in Fig. 7, the momentum equation for the air-water 


mixture can be written as: 


dv dp 
Waeg ha ET eo (6274) 
dz daz 
dp P W dv 

Or Me oo eS (28) 
dz A A daz 


where W is the mixed mean mass flow rate of air and water, V 
is the mean velocity of the mixture within the element, Z is 
the direction of flow, P is the perimeter of the element in 
contact with the pipe wall, A is the cross-sectional area, 
ee is the average wall shear stress, p is the average 
pressure at a pipe cross-section, Pm is the average density 


of the mixture, and g is the gravitational acceleration. 


The three terms on the right hand side of equation 28 


can be rewritten as: 


dp d d a 
oe [ . =| i =| (29) 
az azir lazia lazie 
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Fig. 7 - Momentum considerations of a cylindrical section 


in the vertical drop pipe 
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If J represents the volumetric flux which is the 
volumetric discharge divided by the cross-sectional area, Q 
is the volumetric discharge, V is the velocity, a refers to 
the volumetric fraction of air in the pipe, subscript w 
refers to water, subscript g refers to air, and subscript j 


refers to the total flux, the following expressions can _ be 


written: 
Bocady thd, (30) 
Ya sndg/* (31) 
V, = I/(se) (32) 
een ONE (33) 
Jg = Qg/A (34) 
Vaenny ws he (25) 
Vgi = Vg - J (36) 
Jgj = a( Vg - J) (37) 
oe = nde) - a ) CO Vw. sey) (38) 


where »V yjs.aSeythe drift velocity of water relative to the 
volumetric average, Vg; is the drift velocity of air 
relative to the volumetric average, Jg; is the drift flux of 
air relative to a surface moving at the average velocity, 


and similarly for Jy,. From (37), 


or aj = Jg = aV gj (40) 
Therefore ah=ade/ (4d)* igi ) (41) 
or mesg Ad Fog * Vg ) (42) 


Going back to (29), 
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where M, is the air mass discharge rate per unit area, Mu is 


g 


the fluid mass flow rate per unit area. Hence, 


M 
g 


M 


Substituting (42) 


Vg 


Substituting (42) 


Since 


Vw 


water 


constant down the 


(47), 


dz 
dV, 


dz 


ae 


= PgQg/A (44) 
Py QA /A (45) 
ante: A3rl)s, 
SO ited gy te igs (46) 
imtor se). 
5 Sal Jwtugt Va) ) as 
( JVtVg; ) 


1S “incompressible compared to air, J, is 


shaft, hence, if we differentiate (46) and 


ye S: (48) 


ag J 
wo _ —2/ at (49) 
az 


Also assuming isothermal expansion, 


vg 


Differentiating, 


= ( aaibee /P WS) ee (50) 
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Substituting: (51) “into (52), 


dP 


dz 


J P aP 
4 Yfke | fg ca |: ) fecatm) or! (53) 
) 


g 
A ( IytVg; 


For the frictional component, 


dP 


daz 


2C MJ 
= (54) 
F D 


where Cy DS") the triction factor, M is ithe total mass flow 


rate per unit area, and J is the total volumetric flux, and 


D is the pipe diameter. For the gravity component, 


-dP/dZ), = GPm Coon 
But Bras, Dy at ( 1-a )P. (56) 
Hence, -(dP/dZ), = g( aPgt(1-a)p et (57) 


Subscituting (93).. (54), (57 )into (29), we gets 


dP 2C,MJ/D + glap,+(1-a)oy] 


ens (58) 
GZ pet Ma dy gis ico 


eT 


atm Patm 


This is the one-dimensional momentum equation for the 


two-phase flow in a vertical pipe. 


2.2.2 Comparison with Anderson et. al.'s Data 

Anderson, Vaidyaraman, and Chu(1971) had done some 
observations on the air-water flow problem in a vertical 
drop conduits. They found that the injection of air into the 
pipe is extremely effective in suppressing cavitation. In 
order to compare the analytical derivation with the 


published data of Anderson et. al., equation 58 is used to 
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predict the variation of pressure and velocity. These are 
compared with the actual measurements as shown in table 1. 
The pressure heads versus elevations are plotted on Fig. 8 
for visual comparison. 

The simplified homogeneous two-phase flow theory agrees 
well with the experimental data presented by Anderson, 
Vaidyaraman, and Chu. The application of the theory to 
prototype drop shaft is very limited due to unknown boundary 
conditions at the inlet and the assumption of homogeneity. 
In prototype drop shafts, free surface aeration will 
probably be the dominant inlet condition, which is again 
dependent on the inlet geometry. The outlet pressure will be 
at or above atmospheric depending on backwater condition. 

The two-phase analysis is valid for the experimental 
Situation by Anderson et. al. because a valved inlet 
controls the head loss and is thus closed to the atmosphere. 
The head chamber level is substantially above the inlet so 
that the drop pipe is forced to flow full. Air is allowed to 
enter the drop pipe through air vents evenly distributed at 
one elevation. Thus the rate of air inflow can be controlled 
as related to water discharge and becomes an independent 
variable. However, most prototype drop shafts have open 
inlets and a relatively small and shallow chambers. Hence, 
air flow rate depends on water discharge as well as inlet 
geometry. Flow types become considerably more complicated 
and the homogeneous air-water flow condition is very 


unlikely. 
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Table 1 


Pas) 
Worked example 


Z Az P, Jg Jw J Vg Vw P2 
(kPa) | (m/s) | (m/s) | (m/s) } (m/s) | (m/s) | (kPa) 
1072311453912 .79414-.333 14.803) 4484 | 104.8 
95.6) 15631125794 (4.9424 14.174)45585 1103202 
88.6! 1.75912.794 14255214. 302/4:725) 88.5 


Sie Gps 005 
A OO hee OS 
68. 


Oot 
O39: 
0.409 


P; is the calculated pressure 
P,2 is the measured pressure 


These pressures are plotted in Fig. 8 
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Nevertheless, this analysis may be applicable to cases 
where the drop shaft has been substantially undersized, or 
the discharge being substantially increased by relief sewers 
connected to the inlet. The head chamber is deep enough so 
that during heavy rain storm, the inlet becomes totally 
submerged and the drop pipe flows full. Moreover, the drop 
is so long that cavitation at the neck of the shafts becomes 
imminent. And suppose that some engineer realizes the 
Situation and decides to construct air vents to relieve the 
negative pressure. Under these circumstances, the two-phase 
flow theory can be used to obtain adequate results. In any 
case, it is one of the many possible flow controls for drop 
shafts and its inclusion in the present discussion is 


necessary for a complete description. 


2.3 Flow in Drop Shaft 
If the flow does not fy the entire shaft 
cross-section, it is generally analysed by the Bernoulli 
criteria, assuming the flow remains coherent along the pipe 
wall, and the losses may be lumped into H,,. Hence, 
Vine Vin 


ae Aen 
2g 2g 


ee eee (59) 


where. 6 h Mois: vathe eilengthierdfin thee!dropshaft:, jand \V2 ivs/tthe 
velocity at the outlet of the dropshaft. Therefore, if the 
head losses in the dropshaft are estimated, the velocity in 


front of the impingement area may be computed. 


-_ 
- 

. ee P 

2S2ho O82 

a yi Ssl 


« 
Ai 4 
, 2 ¢ 
i 
rh ® 
' is 
2 ve 
- -_ 
ons «aveii- 
& he e<«a 
Bf: ~OwW 
2 
¥ 4 Mui 
fun « 
Lg « 
mai on 


; : ; . q : ; 
fox Yd heese toni y Let anesedite cried + abuse y 


Oe) aie a sete 
va 7 aL De 
a 7 


mt ah ee 7 a 


‘at ‘sl 7 hi i 
ggeh 2i tedrat> Seep tefnt) ed: ox below 
| a Cs mt. > 37 Oe 
sot dein: sen ieee 
th ; ie 7 fn 
P ype pe baie fr ms boys 
? Sen SB eee Kast #5: tad: pool 
hen we ae i é 
i.208 amOs- Mail) e2emgaue onc or ‘in 
> ee 7 7 wy 
2 v i668, JoompeRoo et 2hhio98 . ng no ai 
s : ‘ys 7 r ene 1 
+ .2ooreyenmesiS BESHE 2eERU yasveseng ovis 
ms | ; a 7 v a y, ey ; 
- an a) a oe a 
. . & Ad 674 + Bee 2c. ngs 4 yoo? : wold 
“ p : un a =* 
- sud : cy 
LOT Koi k SieRE Pee. ae ad 46 | 2am 
: yi , ao ® D ' ay. 
; S229: gts 4 toni ast One, @: $34 
5 . + ore Le 
OLS ceindgHos | & 704 eraeeee 
7 
7 ae 


ey 


"nob toeenagen 


i “ fai B.. = 


28 


2.4 Plunge Pool Hydraulics 

The plunge pool should be designed to effect energy 
dissipation and deaeration of the falling nappe, so that the 
violent turbulence can be confined to the pool area. The 
discharge into the outlet should be quiescent, clear, with 
low air content, and at an acceptable velocity. 

The pool may employ either one of the two energy 
dissipation mechanisms: 1) hydraulic jump; or 2) jet 
diffusion. The geometry of the pool will be distinctly 


different depending on the dissipation principle adopted. 


2.4.1 Hydraulic Jump Basin 

If the pool is designed as a hydraulic jump dissipator, 
it should be elongated with a raised roof similar to the 
pool geometry employed by the TARP project of Chicago. The 
geometry of the energy dissipator is shown in Fig. 9. Weir 
and baffles were not used in the basin because of the 
observed poor hydraulic performance in the model, and _ the 
possible maintenance problems which may arise after every 
severe flood. In Anderson and Dahlin(1975)'s study, the 
basin was found to be effective in releasing the entrained 
air. It also performed satisfactorily for several tailwater 
levels. However, at the design discharge with low tailwater, 
Ehegetiows ins ithe ‘outlet sewer was observed to be 
Supercritical and the estimated outlet velocity was 
approximately 7 m/sec. This velocity is probably too high 


for debris-carrying wastewater. Therefore, a deep sunken 
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Figure 9 - Drop structure geometry adopted by TARP (From 
Anderson and Dahlin(1975)) 
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basin may be needed for more effective energy dissipation at 
low tailwater depths. A forced jump is preferred to a free 
jump because the former requires a much shorter length and 
thus a lesser excavation volume. 

A suggested pool geometry is shown in Fig. 10. As the 
downward flow strikes the plunge pool, it will be deflected 
into a backward flow anda forward flow if a backpool is 
provided. This has the desirable feature of strong 
circulation behind the flow which assists in energy 
dvssipatiom.) As, a _.result,. the kinetic energy of the 
Supercritical flow in, front of the jump may-well be less 
than that predicted by the potential flow theory and the 
actual required basin length can be shortened. The benefit 
of providing the backpool and the potential for cost saving 
in a design Situation will of course require the 
verification by a physical model study. 

Upstream of the outlet sewer, a ground roller will form 
where debris will be trapped in circulation. This has_ the 
harmful potential of scouring the basin floor. A possible 
solution is to provide a sloping floor at the sewer entrance 
so that the basin will be self-cleaning and the erosion 
problem minimized. 

Considering the control volume A in.Fig. 8, the force 
balance in the horizontal direction may be written as: 

F, = Fs + Ms; (60) 
where Fy is othe. back pressure, force, 9F3 1s» the total 


pressure at the y; section, and M; is the momentum flux 
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Figure 10- A suggested geometry for the forced jump basin: 
(a)sectional view; (b) plan view 
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exiting the control volume A. 

The forced hydraulic jump in a uniform rectangular 
sunken basin was analysed by Rajaratnam(1964). The effect of 
the step was accounted for by introducing a drag force, P,, 
into the momentum equation. The drag force was defined as 
[Cu (pV3?/2)-Ah], where Cy is the drag coefficient, and Ah 
is the basin step height. With this, the momentum equation 


for the forced hydraulic jump may be written as, 


Yrs YYo’ Yq 
Beeee ir: =—-*[ Vor - V3 ] (61) 
2 2 g 


where y3 is the supercritical flow depth in front of the 
jump, Yo is the flow depth in the outlet sewer, q is the 
discharge per unit width, Vo; is the outlet velocity, and V3; 
is the supercritical velocity in front of the jump. 

In Rajaratnam's experiments, at least six distinct 
forms of forced jumps could be observed by varying the 
distances off theystep -from the, stant 4.0fsi thesagump,4 Xon.\y in 
order to minimize excavation volume, the basin length should 
be kept as short as possible. Nevertheless, energy 
dissipation and jump containment must remain effective. The 
values of Cy were determined experimentally for various 
Xo/L,; ratios. A design chart was also developed so that for 
a given supercritical flow with the known tailwater level, a 
Suitable basin length and step height may be chosen. 

Presumably a similar procedure may be used to design 
the pool in Fig. 8. Some experimental observations will be 


needed to develop the design criteria and finalize the 
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appropriate pool geometry. 


2.4.2 The Jet Diffusion Mechanism 

If the pool is designed to effect energy dissipation by 
jet diffusion, it is recommended that its shape be circular. 
The shaft is preferably located directly above the center of 
the plunge pool so that the jet may expand in all 
directions. A sketch of the possible flow patterns in the 
pool is shown in Fig. 11. Due to the presence of the pool 
boundaries, the flow in the plunge pool may be divided into 
three regions: 1) the circular downward jet region; 2) the 
impingement region; and 3) the annular upward wall jet 
region. The extent of each of these regions will depend on 
the initial momentum of the jet , the jet size versus’ the 
plunge pool dimensions, air content, and pressure 
distribution in the pool. 

The diffusion mechanism of a deeply submerged circular 
jet was discussed by Rajaratnam(1976). Since the jet belongs 
to the class of slender flows, it may require a pool depth 
of at least one order of magnitude larger than the jet size 
in order to achieve the required level of energy 
dissipation. As a result, it may not be economical to design 
an underground stilling pool utilizing the jet diffusion 
mechanism, unless some submergence exists in the _ plunge 
pool. Besides, debris is likely to be collected at the 
bottom of the plunge pool which will require routine 


dewatering and cleaning. 
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Plunge pool 


projected area 
of Outlet sewer 


downward flow 


Downward flow 


Circular 
jet 
region 


depth 


Impingement region 


Figure 11- Possible flow patterns in plunge pool: 

a. when the pool depth is deep enough and the 
pool diameter large enough to effect energy 
dissipation 

b. Plan view of plunge pool 

c. when the pool depth is not deep enough to 
dissipate the kinetic energy of the flow, and the 
jet is deflected into the outlet sewer 
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Without any experimental observations, it is difficult 
to determine the most efficient plunge pool geometry. For 
the purpose of comparison, the plunge pool previously 
adopted by the City of Edmonton is’ evaluated using the 
forced jump criteria. This plunge pool geometry is shown in 
Fig. 12. A typical shaft size of fifteen feet is selected 


for the discussion in the following section. 


2.4.3 Evaluation of the Plunge Pool Geometry Using the 

Forced-jump Criteria 

To evaluate the adequacy of the plunge pool, a_ typical 
shaft drop structure is chosen for analysis using the forced 
jump criteria. Since the pool geometry has not been tested 
experimentally, there is no data available which would 
indicate the possible flow characteristics and the proper 
length requirements. In order to proceed with the analysis, 
it 1s assumed that the circular basin can be replaced by a 
Square pool with the same major dimensions. 

tia, pager cular ‘drop shaft has the following 
dimensions: 


D Putts 


h 1.0 Ot. 


then, the plunge pool size is 1.5D or 22.5 ft., and the pool 
height. Ah, is. 8.22 ft... Tf the inlet sewer is horizontal 
and is designed to have a y,/D value of 0.5 at the design 


discharge, using (14), the design discharge is computed to 


be 1485 cfs. 
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Figure 12 - Plunge pool geometry adopted by the City of 
Edmonton 
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A wide range of tailwater levels may exist depending on 
the downstream conditions. The most critical condition is 
probably when the outlet sewer is unsubmerged. In this 
Situation the minimum tailwater depth, yo;, that would occur 
is probably equal to the critical depth. If both the inlet 
and the outlet are of the same size, this flow depth is 
equal oto ./35 ft.\., From. (26), 
AEN ito h Ets, 

from (24), Vera ts 5. Matt. /Sec.. 

and from (59), neglecting the losses in the shaft: 
V5 = 63.3 .ft. /Asec. 

Now, 1f~the velocity. in «front ~of the jump, V:, is 
assumed to be the same as V2, and the pool width, B, is 
equal to 1.5D, the average flow depth in front of the jump 
is 0.81 ft., the supercritical Froude number, F3;, is 16, and 
the sequent flow depth, yz, which would exist in a 
Norizontvalsoasini1s 174.9 st... Hence, 

Vaz Vo/ Va =H) o/ 17.9 pe 0.42 
The distance from the step to the toe of the forced jump, 
Ke, iS) sapproximatesy 0.75D. ~The jlength of the »surface 
roller, Lei , vs approximately 4.5 times the sequent flow 
Gepth. Hence, 

Ko/L,; = 0.14 

Brom Figs 13, °C. is 0.65 andi the “gumpiyeis within the 
type IV region but rather close to the type VI region. Type 
VI jump is not acceptable because the supercritical flow 


simply shoots over the step to forma swell behind the 
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basin. This indicates that the basin is probably too short. 

From Rajaratnam(1964), 

( PY aI OAD ReAeE ow. f LOW ye 
F37ov, 

where A = BCug, @ = y2/Ys, Ya = Yoi¥ay B ="Ah/y.3, Ca is the 


(62) 


drag coefficient due to the basin--step, jy, is the 
Subscritical sequent depth, y; is the supercritical flow 
depth in the basin, and yo is the outlet depth. For the 
present drop shaft, 

A = BCa = 1.45 
Hence, Be= Dnfyer= 2.23%, 
and An (required) = 1.81 ft. 

Therefore, the Ah(provided) is approximately 4.5 times 
the required pool depth. If the basin is sufficiently long, 
energy dissipation will remain effective with a submerged 
hydraulic jump formed in the basin. However, if the length 
of the basin is too short, the supercritical flow will 
Simply be deflected to enter the outlet sewer. These 
possible flow patterns are shown in Fig. 14. 

Rajaratnam(1964) suggested that for good hydraulic 
performance, a type II to type IV jump should be chosen. For 
the following discussion, a type III jump with a drag 
coefficient of 0.4 is selected. From Fig. 11, 

Xo/L,j = 0)..53 
Hence, Xou=442, 8 ft. 
and the minimum pool diameter required is: 
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surface roller 


ground roller 


Supercritical 
flow discharged 
into 

outlet 
sewer 


Figure 14 - Possible flow pattern in the hydraulic jump 


basin: 

a. when the step is bigger than that required to 
form a forced jump but the length of the basin 
is long enough to contain the submerged jump; 

b. when the step is deep and the length too 
short to form a jump, the supercritical flow is 
deflected into the outlet sewer. 
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A backpool may be provided if model studies indicate the 
benefits of incorporating it into the structure. 

Since A is 1.451, and B is 3.63, the required step height is 
2.9 feet. 

The above calculation provides an approximate method of 
evaluation for the plunge pool geometry shown in Fig. 10. 
However, for future design, a model study is required to 
determine the adequate geometry and to verify the hydraulic 
concepts discussed. Hydraulics of plunge pool can be 
modelled fairly accurately using the Froudian criterion. In 
most cases, substantial reductions in construction and 
maintenance costs can be achieved by a laboratory model 


study. 


2.5 Effect of Outlet Submergence 

When the tailwater level increases above the sewer 
crown, the outlet sewer may flow full. In this’ Situation, 
the water depth in the shaft may be determined by uSing the 
Ppipe-flow equation; 


L H A Vg 
Hy = [ie hog et frat fy'—: —] |: (63) 
D Dy 2g 


where Kg is the entrance loss coefficient, Kp 1s the exit 
toss laccett tement Rik, vers. ithe Darcy-Weisbach friction factor 
for the dropshaft, L is the length of the outlet sewer, H, 

is the submerged length of the dropshaft, D, is the diameter 


of the dropshaft, A is the area of the outlet sewer, A, is 


the area of the dropshaft, and Vg is the full flow velocity 
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of the outlet. 
itsenefouclet conduit is horizontal, thé depth in ‘the 
vertical shaft, H;, can be calculated as: 
Bape title (Ney com) (64) 
Orit ythe outlet conduit is at a slope 'S,, 

Deere He Sob: =. (7 O2/2qhezes) (65) 
where T, is the tailwater depth, and Mee is the average 
velocity in the drop shaft. A sketch of the submerged plunge 
pool is shown in Fig. 15. Submergence of the outlet will not 
alter the discharge capacity of the drop shaft unless’ the 
water level in the shaft rises to drown out the critical 


flow control section at the inlet. 


2.6 Overloading 

Since the discharge capacity is controlled by the 
terminal weir at the inlet, it will not be altered unless 
y,/D becomes greater than 0.9 or the shaft entrance is 
flooded out. In the former case, the inlet conduit will flow 
full except near the outlet. The driving head is the head 
difference between the upstream manhole and the shaft 
entrance, minus the friction losses. The head discharge 
relation may be described by the full pipe-flow equation. In 
the latter case, the vertical shaft is completely filled. 
The discharge capacity may be determined by the full flow 
capacity of the outlet. The discharge in the inlet should be 
the same as that of the outlet under steady state 


conditions. If the inlet discharge thus computed is higher, 
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the water level in the shaft will rise to increase the 
discharge in the outlet. Otherwise, the water level will 


fall and the outlet discharge diminished. 


2.7 Uplift Pressure 

if Mthel drop structure is Situated below the 
groundwater-table, the pore water pressure at the base of 
the plunge pool may be sufficient to heave or crack it. The 
most critical condition occurs when the structure is empty. 
If field tests indicate that seepage forces will be a 
problem at the site, control measures such as providing a 
gravel blanket with drain pipes beneath the basin floor 


should be considered. 


2.8 Air Entrainment 


2.8.1 Types of Entrainment Processes 

There are possibly four different kinds of air 
entraining processes which can be identified in a vertical 
drop shaft. These are: 1) the air entrainment due to. the 
Surface roller at the shaft entrance; 2) the free surface 
aeration of the flow in the shaft; 3) the forced aeration as 
the flow enters the water surface in the shaft; and 4) the 
air entrainment due to the hydraulic jump in the stilling 


pool at low tailwater. An assessment of these entrainment 


processes are as follows. 
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2.8.1.1 Surface Roller 

When the shaft entrance does not have a smooth 
transition, the nappe will leave the inlet invert to impinge 
on the outer shaft wall. The flow will be deflected downward 
with a surface roller formed on top of the upper nappe 
causing local aeration. This aeration phenomenon was studied 
by Renner and had been summarized by Rao and Kobus(1975). 
Although the surface of impingement in Renner's study was 
not curved, the flow geometry is generally similar to the 
present case. In the absence of model studies, his results 
could be used for a rough estimate. 

Some of the results in Renner's experiments are shown 
in Fig. 16. He showed that if the concentration of air in 
the deflected nappe beneath the surface roller is far from 
the state of saturation’, the rate of air discharge to water 
discharge is proportional to the square of the jet's Froude 


number®?. 


2.8.1.2 Surface Aeration in the Drop Shaft 

As the flow runs down the shaft, its free surface will 
entrain air causing the flow to expand and become foamy. If 
the shaft is long enough, the jet may disintegrate into 


droplets and slugs of water. From limited experiments* 
2 The state of saturation corresponds to an air 
concentration, C, of about “0.4 to’ 0.5, Or an air to water 


discharge ratio of about 0.7 to 1.0. 
> For aeration to take place, the jet's Froude number must 


be greater than 1. 
* The range of their experiments 1S: 
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Figure 16 - Relative air entrainment in a horizontal jet 
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McKeogh and Elsaway (1980) established the following 


relationships for circular free waterfall: 


For ewer 0.80%,0L, = 66> Oar” "4G (66) 
£Or e = 1%, Ly = 3820 2a. ae (67) 
and for ceet 5%, hyp = 436: OURPEa, (68) 


where Ly, is the disintegration length, Opnds the water 
discharge, e« is the turbulent intensity defined as (/u'?)/u, 
and u, is the jet velocity. 

The process of free surface aeration in a _ rectangular 
flume was investigated by Straub and Anderson(1958). They 
Studied the air concentration due to self-aeration in a 
rectangular flume with the slope ranging from 7.5° to 75°. 
From their experiments, they found that the mean air 
concentration in the fully-developed flow is a function of 


the slope and water discharge: 


£4. (Shq 0 (69) 


Cm 
where Cp, iS the mean air concentration of the fully 
developed flow, and S is the flume slope. However, it is not 
known whether their results can be extrapolated to represent 
the vertically flow in saodropo%shaft,’Itcis likely>thatrthe 
process of aeration is similar in the early stage of the 
fall, with the characteristics deviating significantly as 
the jet begins to disintergrate. Again, further research is 


needed to clarify this. 


*(cont'd) Os7O pie tet 0 > m/sec. 


Note that the disintegration length, Ls, is a strong 
function of the jet's turbulence level. 
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2.8.1.3 Forced Aeration by a Plunging Water Jet 
The mechanisms of aeration when a_ high velocity flow 

enters a slower moving body of water was classified by 

McKeogh and Elsaway(1980) as: 1) annular oscillation; 

2) intermittent vortex; 3) turbulent occlusion; and 

4) droplet entrainment. Detailed descriptions of these 

mechanisms may be found in their paper. 

AS. the swater »ilevel..ine» the «sShafit. rises,.¢ther air 
transport capacity may vary. The rate of air transport will 
be governed by the entrainment capacity of the falling jet 
as well as the transport capacity of the flow below the 
water surface. Willock and Thorn(1973) identified three main 
regimes of air demand: 

Ave When the nominal fluid velocity is less than the 
threshold velocity for air flow, the entrained air will 
be retained in the zone of recirculation. Only very 
fine air bubbles can be discharged into the outlet 
tunnel below. Hence, the air demand is almost zero. 
This threshold velocity was determined to be about 0.13 
m/sec. which is about one-half of the bubble rise 
velocity, for the range of bubble sizes between 1 mm 
and 5 mm in diameter. 

2. In this case, the nominal fluid velocity remains less 
than the threshold velocity but the zone of 
recirculation is not completely contained within the 
shaft. As a result, some air is carried into the outlet 


tunnel. The length of the recirculation zone was found 
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to be about five to six times the shaft diameter and 

air transport will occur if the water column in the 

shaft is less than this value. 

a When the nominal shaft velocity is greater than 0.5 
m/sec., which is about twice the bubble rise 
velocity®, all the entrained air will be transported 
downward. The water level in the shaft will not change 
the air demand characteristics and the ratio of air to 
water discharge becomes a unique relation of the fall 
height. The air discharge is caused by a combination of 
self-aeration of the falling jet and the entrainment 
during impact with the water surface. The experimental 
results of the discharge ratio versus drop length is 
shown in Fig. 17. 

In Thorn and Whillock's experiments, they found that 
the discharge ratio reaches a maximum at fall height of 
about 3.7 m. However, McKeogh and Elsaway(1980) found that 
the maximum volume of air entrainment is consistently about 
10% higher than the volume at a fall height equal to the 
disintegration length. Since the disintegration length 
depends on the turbulence level and the fluid discharge of 
the flow, there is an apparent contradiction between the two 


investigations. 


> The majority of air bubble sizes range between 1 to 5 mm 
in diameter. In this range, the air bubble rise velocity is 
about 0.26 m/sec. This size range has been found to be true 
for both models and prototype structures. 
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2.8.1.4 Air Entrainment by a Hydraulic Jump 

When the tailwater level is low, the jet will strike 
the plunge pool at the bottom of the shaft where further 
entrainment takes place. If the pool is long enough to allow 
the formation of a stable hydraulic jump, the entrained air 
will be released by the surface roller. In this case, air 
demand will not be a problem as long as the plunge pool is 
ventilated. The rate of aeration by a rectangular hydraulic 
jump is a function of the Froude number which was shown by 
Kalinske and Robertson(1943) to be: 

Be =505,0066-1. Fy=4 ie) 7g (70) 
where B is the ratio of air discharge to water discharge, 
and F is the supercritical Froude number. However, if the 
pool is.tooysmall.to,.contain,. the ,jump and to allow the 
entrained air to escape, it will be discharged into the 


outlet tunnel. 


2.8.2 Air Entrainment Modelling 

Due to the complexity of the entrainment phenomenon, 
there is no satisfactory theoretical treatment available at 
this time. For the design of large hydraulic structures, 
models are usually made to investigate their operational 
characteristics and the various design options. Most dynamic 
aspects are modelled satisfactorily according to the 
Froudian criterion. However, such model studies tend to 
underestimate the rate of air entrainment. Kenn(1967) 


Suggested that the inherent scale effect in predicting air 
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entrainment using the Froudian model may be due to. surface 
tension and viscous effects. Ervine(1976), and McKeogh and 
Elsaway(1980), also added that in the case of air 
entrainment due to a water jet entering a pool of water, the 
turbulent intensity of the jet bears an important effect on 
the amount of air entrained. 

Kenn(1967) proposed that if the model velocity and 
pressure force are of the same order of magnitude as that in 
the prototype, air entrainment may be Satisfactorily 


modelled. Considering viscous and surface tension forces: 


Weber number pV?1/a uv 
ee C71) 
Reynolds number pVl/u 0 


which requires that (uV/o) is the same in both the model and 
the prototype structure. Using the same fluid in the model 
as in the prototype requires that V(model) be the same as 
V(prototype). 

In Ervine's(1976) study, he found that for a jet 
Striking an ambient water body, the ratio of air to water 


discharge is given by: 


a 


fe) b H O.446 Vio 
Boa! © |-[ sian (72) 
Ors P a V 


where K is the proportion of the nappe perimeter exposed to 
the atmosphere, b is the nappe width, d is the nappe 
thickness, V is the velocity of the nappe, Vp is the 
incipient entrainment velocity, H is the height of fall, and 


P is the nappe perimeter. 
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The equation also indicates that air entrainment cannot 
be modelled by the Froudian criterion but requires that 
(Vo/V) be the same in both the prototype and the model. As 
indicated by Ervine, the incipient entrainment velocity is 
approximately the same for various jet sizes. In order to 
maintain similarity between the model and the prototype, the 
nappe velocity must be kept the same, or else the (V,/V) 
ratio must be very small compared to one. Whillock and 
Thorn(1973) also concluded from their experimental studies 
that a distorted model with increased drop length, to allow 
the downward flow to attain prototype velocity, may provide 
a reasonable prediction of air demand. 

The knowledge of air entrainment in a drop shaft is 
necessary in the design of the ventilation system. Previous 
experiences have shown that serious operational difficulties 
and a loss of service may result if air entrainment is not 
accounted for properly. Unless extensive experimental work 
is to be performed, the air demand of such Structures may 


best be estimated from field measurements. 


2.9 Modifications 


Based on the above discussion, the following 
modifications may improve the efficiency of the shaft drop 


structures: 
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2.9.1 Inlet Transition 

The entrance to the drop shaft should be _ smoothly 
curved. The flow cannot change its direction instantly and 
the presence of abrupt directional changes will cause the 
nappe to leave the lower boundary. The formation of a 
surface roller by the flow impinging on the pipe wall was 
discussed in an earlier section. This is undesirable because 
of the enhancement of turbulence and the reduction of the 
dropshaft capacity. 

Flow around a vertical elbow may be roughly explained 
as follows: 


Referring to Fig. 18 and assuming irrotational flow, 


Kile vr = Vo LOG ya eTeifAs se G73) 
Bony 
and Ov= \ © (vB)dr (74) 
Yo 
where B = Dsiné, (75) 
and dr = 0.5Dsinéde (76) 


Subsitute (73) into (74), we have: 
cos” “Gi=2y,/D) 
K/2 
[——___|p'sinve dé (77) 
rot+tD(1-cosé)/2 

0 

For a given discharge and curve radius, K can be solved 
for by using (77) and the velocity at any depth computed. 
Due to the increase in velocity near the boundary, the 
pressure may approach a cavitating level at high discharges. 


Therefore, the inlet transition curve should be smooth with 
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Section A-A 
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Figure 18 - Smooth inlet transition 


‘+ " 7 
, i i ; y ea 7 
i] 
=) 34i o*@ 7 : i i" ; 
é . 
: a5 Lm | 
: 7 ‘ Say 
a) 
a rn i 
wl i o 7 
ae ss 7 
™ ; 1e » j F 
i | oe 
a * , 
—— er oe OD i 2 -_ — 
- ’ ~~» 
bd e : 
Vy ¢ — a s ate | 
a ‘ 
eee | 
' i i ‘ 
- a | 
¥ : ‘ 
: « nn a — | : 
- * 
. j i 
4% i 
: 
| ~ ! 
i ~> o 7 7 


; | “ase, ren 
a peal | in ee as Ly 
4 ' 


a 
\ rT a ee F Na 
a e’ a lai. 


ay =-7 


56 


a gradual directional change, or else the lower nappe should 


be ventilated. 


2.9.2 Shaft Size 

The drop shaft size should be at least the same as’ the 
inlet sewer size. Since there is no accurate way of 
predicting the air discharge, the higher shaft capacity can 
be used as a safety factor against unforseen conditions. In 
Anderson and Dahlin(1975)'s experiment, performance of the 
drop shaft does not seem to be impeded by the insertion of a 
divider wall other than a mild increase in the impact 
pressure measured in the plunge pool directly below the drop 


Shaft. 


2.9.3 Divider Wall 

The partition wall for the inlet air vent should be 
slotted to provide continual aeration of the downward flow. 
As the flow accelerates down the shaft, its air demand will 
increase continuously. Tt air passage is impeded, 
subatmospheric pressure may develop. This suction head can 
cause a violent insufflation of air through the entrance. 
Besides, local pressures may reach the cavitational level 


cauSing excessive erosion to the pipe surface. 


2.9.4 Plunge Pool Sizing 
The possible modes of energy dissipation mechanisms and 


the appropriate plunge pool geometries were discussed 
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previously. Generally, it is desirable to have an elevated 
roof in the pool to allow the passage of returning air into 
the inlet air vent. A longitudinal pool is preferred over a 
circular pool in the TARP experience, which is also equipped 
with an inclined roof to provide ample clearance for air 
discharge. Venting the discharged air through the inlet will 
have the advantage of recirculating the air in the system 
and hence reducing the air demand from the inlet chamber. 

If the pool is not long enough to contain a well formed 
jump, the supercritical flow may discharge into the outlet 
tunnel. This not only causes excessive erosion problems’ in 
the basin and the outlet sewer, but also allows the 
entrained air to be discharged into the outlet. In this 
Situation, a downstream air-vent will be needed to prevent 
the formation of high pressure air pockets along the roof of 
the outlet. These air pockets may blow back through the 
entrance and set up unsteady surges in the shaft. 

The optimal plunge pool geometry can best be determined 
from a model study. A forced jump basin with the features 
illustrated in Fig. 8 would probably be adequate for both 


energy dissipation and deaeration. 


2.9.5 Outlet Sewer Sizing 

From Anderson and Dahlin(1975), an outlet size equal to 
the inlet size is adequate to handle the design flood if the 
impact basin is well designed and ventilated. Indeed, if the 


outlet size is intentionally reduced, the flow depth at the 
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design discharge with no tailwater will be increased. Hence, 
the plunge pool depth may be proportionally reduced, 
resulting in a more economical basin design. However, when 
the tailwater level rises, and the shaft and the outlet 
conduit flow full, the energy losses from the entrance to 
the outlet may become critical. A smaller outlet size will 
have higher energy losses and a lower discharge capacity. In 
this situation, the entrance to the conduit may be bevelled 
to improve performance and to reduce head loss. In order to 
achieve optimal sizing of the outlet conduit, the range of 
Operating tailwater must be known, and the plunge pool must 


be sized to effect energy dissipation and flow deaeration. 


2.10 Recommendations 
Based on the previous discussion, the following points 

can be suggested for the design of shaft drop structures: 

Pe The flow characteristics in a smooth entrance is 
Superior to an entrance with abrupt changes. Various 
entrance configuratons should be tested by a model 
study to obtain an optimal design between economy and 
efficiency. 

25 As discussed earlier, the flow conditions in the drop 
shaft can be very complicated and are heavily dependent 
on the entrance geometry. The actual flow 
configurations, be it a self-aerated open channel flow, 
Or an air-controlled flow, and the possibility of local 


Cavitation and separation in the shaft can best be 
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studied using a transparent model. 

The theoretical derivations and assumptions presented 
should be verified by experiments’ before their 
application in future use. 

A suggested plunge pool geometry is provided 
previously. The adequacy of the plunge pool should be 
tested with a model. 

The outlet sewer does not need to be the same size as 
the drop shaft. The suitable size may be determined if 
the range of operating tailwater is known. 

Extensive research has indicated that air entrainment 
cannot be modelled by the Froudian criterion. Some 
field measurements on the operating structures will be 
useful in sizing air vents which may be used to 


Supplement experimental observations 
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3. Drill Drop Structures 


3.1 Hydraulics 

The hydraulics of the drill drop structures will be 
discussed in two parts: 1) the inlet sewer capacity; and 2) 
the drop Pipe capacity. The efficient operation and 
structural economy will depend on the balanced sizing of the 


two components. 


3.1.1 Inlet Sewer Capacity 

Inlet Sizing has been discussed in chapter 2. 
Generally, the inlet should be designed to operate with free 
Outlet condition. lf the conduit slope is mild, its capacity 


may be determined from the terminal-weir-control equation. 


3.1.2 Drop Shaft Design 

To determine the required sizes of the drill drop pipe 
and the head chamber, the entrance capacity must be known. 
In other words, the head-discharge relation of the 


horizontal orifice has to be determined. 


3.1.2.1 Previous Studies on Horizontal Circular Weirs 
Extensive investigations have been done on _ the 
discharge characteristics of circular weirs. Discharge 
coefficients had been given by Wagner(1954), Camp and 
Howe(1939), and Gourley(1911) for various weir heights and 


diameters. 
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Binnie(1938) had made some careful observations on the 
flow over a re-entrant vertical tube. From his studies, he 
classified the flow over circular weir into four possible 
categories: 1) weir flow; 2) Borda full flow; 3) Borda free 
flow; and 4) full-pipe flow. A description of these is 
summarized in Appendix III. 

Although there is an abundance of research on circular 
weir, research on weir with zero sill height, which is quite 
common in many smaller structures, seems to be lacking. In 
the following section, the discharge over circular weir with 
zero sill height is derived analytically and the result is 


verified by experiments in the laboratory. 


3.1.2.2 Circular Weir with Zero Approach Height - Momentum 
Analysis 

Referring to Fig. 19, the momentum equation between the 
critical section and the end-depth section may be written 
as: 

Bip 2h sin (872) 2K Po pOUAV we) (78) 
where Be is the critical pressure at sector A-A, P, is the 
Side pressure force, K; is the pressure coefficient, Pe is 
the end depth pressure, Vo is the critical velocity, Vomis 
the end-depth velocity, and 6 is one-half of the inscribed 
angle of sector A-A. In the equation, 

6/2 
P = ( p.cos#r. de So ee (79) 
=O /e2 


where g is a dummy variable, 
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and Pp Ve vy,/2 
Integrating with respect to the limits, we get: 

PEPE, yyie xr, sink e/2) (80) 
and similarly, a= vy °R-sin(6/2) (81) 


where r. is the critical radius, and-R-is~the-radius of ‘the 
drill drop pipe. 

For an approach slope of zero, K, may be taken as zero. 
Therefore, the third term in (78) may be neglected. For the 


Side pressure force: 


yy 
P, = | dr (82) 


If the water surface profile is approximated by a straight 


lanes and fer Ris r & 256 


Cc 
Ac aes Cy (oy eR) AGre sR) (83) 
Substituting x= foeR) (ir oR), Gras (oak ax and (83) into 
(82), 
1 
y 
P, = . | Dye ay (Cy oe) yey ye A Gogh) dx (84) 


0 
Evaluating the above integral, we get: 
Prem acy homme rR isl sya tay ven ted (85) 
If the discharge in sector AA is given by: 
6/2 
Qe = fossa y,V,cos#r. de (86) 
The critical velocity can be written as: 


N= Oo/l penny usec). (87) 
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Similarly, the velocity at the end section is: 
VS ae Qo/[ 2Ry ,:sin(6/2) ] (88) 
SUDStIrtutang "sO GB), “C85.) 91087», fand’ (88): into: C78.) 


v¥,°r.7sin(6/2) = Racal ves + y si dasa gi Jsiné 
Peni eae a ince a 
2Ry ,°sin(6/2) 20 V eeoun 9/2) 


e 
(89) 


Equation 89 can be rearranged into: 
2gsin AGO /2): = (2/3)gsin? (0/2)-f1+y./y.+ly,/y,)?7]+(1-R/r,) 

= Oe" /Yaet. “NC /RyS) ACT AS ye (90) 
Now, Pit? we'tsubstrtute f).15° ,| sin 6 /2)8t=00), 06719 10 = 12048; 
and g = 9.81 m/sec? into (90), and assume that (r.-R) = ay 


andi(yey/y,, i t= t0in7 251% weieget : 


: Has Say. Rt4y . 3 944y | REayS 


V 0.5256-( 0.275 + 4y /R ) 


(91) 


Rearranging into a non-dimensional form: 


Ry AATEAS PCA SiS PL IR EAS PLS 
Vv (10.275 yee) 

where Q* = Q/g'/?R5/?2 

Equation (92) gives the relationship between discharge 
and end depth for a circular weir with zero approach depth. 
Experiments were performed in the Graduate Hydraulic 
Laboratory to verify the above derivation. Details of the 
experiments will be discussed in a later section. 

Let the head in the chamber be H: 

Pvt fat ( V,°/2g ) (93) 


Rearranging into a non-dimensional form: 
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y 1 1 
ay -Q*? -—____ (94) 
R 81? (y,/R)? 


Combining (92) and (94) will give a head-discharge relation- 
Ship which is more convenient to use than (92) alone. 


THis relatvonsitp—is plotted in “Fig. ‘20% 


3.1.2.3 Full-pipe Flow and Cavitating Flow 

As the head in the chamber increases, a point will be 
reached when the entrance water surface becomes level and 
the pipe begins to flow full. It is assumed that orifice 
control will not occur and the requirements of such is 
discussed in the next section. Then, depending on the length 
of the drop pipe to its diameter, cavitating flow may or may 
not occur. 

The equation for the full-pipe flow is: 

Hee Gomes fei Lig /D i) tei eon) ena (95) 
where H is the head above the inlet entrance, H, is the head 
loss at the entrance, and Lo is the length of the drill drop 
Pipe. If the entrance loss is neglected, (95) may be 


rearranged into; 


Vi R R 
OR fe (96) 
/ Lo 
‘a 
V D 
The equation for cavitating flow is: 
BGP)! “eae (V?/2g) (97) 


Where pers) is the cavitating pressure, and H, is the lumped 


head loss. Equation 97 may be rewritten as: 
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VE H P Hy 
Or = /2n*| —-- = - — (98) 
V R RR 


For an open system with a head chamber above the drop shaft, 
and the entrance is designed so that orifice flow will not 
occur, three types of flow are possible: 1) the weir flow; 
2) Peete Cavitating flow; (and 8) esthe full-pipe flow. 
Cavitating flow has been described by Anderson, Vaidyaraman, 
and Chu(1971). It occurs when the drop-pipe is long enough 
so that the pressure in the upper part of the conduit 
reaches cavitating pressure, Pou husmitype of flow is 
Particularly undesirable because cavitation will erode the 
shaft surface rapidly and severely. 

Drop pipe can be designed so that cavitating flow will 
not control at any stage of discharges. From (47), the 
discharge for full-pipe flow is given by: 

Omres f.( H/R; Uo/Ri Hey Ry flee) (99) 
and from (98), the discharge for cavitating flow is given 
by: 

Overt; (<H/R: (Dy Hak.) (100) 

Hence, if the roughness of the pipe and the loss at the 
entrance is estimated, cavitating flow can be avoided by 
choosing the length to diameter ratio properly. Equations 96 
and 98 are plotted in Fig. 20 for several sizes and L)/R 
ratios. The size of pipe should be chosen so that the weir 
flow curve will intersect the full-pipe flow curve before it 
intersects the cavitating flow curve. In any case, the drop 


Pipe should be designed to operate within the weir flow 
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regime. Pipe flow would occur only in a state of surcharge. 


3.1.3 Computer Program 

The various’ flow controls discussed above are 
consolidated into a design program to facilitate the design 
of this structures. The program will evaluate, for a certain 
drop pipe size, drop height, and outlet water level, the 
head in the inlet chamber and the corresponding discharge. 
This is convenient LOG the evaluation of existing 
structures. In a design situation, the design discharge and 
drop height is known, the drop pipe size and the head in the 
inlet chamber are to be determined. In this case, one can 
determine the head capacity curves for a range of drop pipe 
sizes. The suitable drop pipe is the smallest size which 
passes the design discharge within the weir flow regime. 

The program listing is included in Appendix III. Sample 
output! ris) salsoiincludedfor ‘a typical run. Due tothe fact 
that head discharge relation deviate from the weir flow 
curve as the discharge increases, allowance should be 
provided to the head chamber depth so that surcharging will 
not occur prematurely. Besides, the program does not 
eliminate the need for a physical model study to verify the 


no-orifice flow assumption. 


3.1.4 Methods to Prevent Orifice Flow 
As presented earlier, orifice flow may be prevented by: 


1) a smooth inlet transition; or 2) an enlarged entrance; or 


iva wig - We = 
2. °. 7 ? 7 a 


sepratosue 36 sfaca 6 Ti. Lao RUT | ‘biwew wel? eget 


i ’ T i 
rT , ‘ i | : 
i ave wal wa a 
7 
aw ap1g03 is 109 ee 
j — | ot a 


=| : t => ae i oa 
= "a8 > ; 7 w : { LO z9nAG3 *) ose ; #ihc ag’ wi 


mete a. e ‘ 3 a4; : .> = : if me TP | ap tae S ¢ hi b 


MipekG $A? '.aorud st 
~ “HP tad | F Pa 
‘ si oe ite? 2a? va ese 
' , 
- 7 z. 1 i ‘i eh 2h irs 
RY a 2 its oo 
4 usia npiaeh a vit 
/ 
; I 1 Silt , wee c 
. : igJen o4. gv =% 
e 
f 1 RSoBgeH Pe sng 
~ ‘ i Ms z *c4g0 wot. fees bea 
a - ye * Fo, ls i) LW ol irhspadelac mie 5 
> a iy * - r q hy a eri * - ro 
a Oneal bid Tous \ bh OTL £ Riana ma Tyem be a 
ee \ Nase 
nel sd404. sud, fo4 Lestaye ele BAbiipans oma at fuq@3t 
' arn a3 ‘ Ty ae 
7 r 4 ane esos s ") ari2 
me ia a ) Bie : ae a~! 4 4 od ob wabeee 3 ‘ag intiow bh Sead 908 
eo 


‘S -W.Vens - songyolle oe eotinchsniti caliph: ra 


Sy ey! 
— ~« tk mee > ae wy 
‘ PI SET 2 740} Ce sae “Stdiins dinette. bite aise 


| i G>, per iatecoh 
: . ) aise m 
wid yiitev'os ybude: behon- =}2 aqaq wheel ia 
oa aay be ts 
citi 


en 
7 ey. | : : pak) 2”. 
sed Se peenaven: ad) Yan Baie 10 \SBets 


7) ur) “- sat . 


T rae ae bap rsine: fp 


7 
a re 


—— 


69 


3) the length to diameter ratio kept larger than 24. 


1. 


The inlet may be designed to follow the lower nappe 
protilew-on Circular ‘overfal) similar to that ofa 
morning glory spillway. This entrance is efficient with 
little loss and no cavitation problems. For small 
Structures, a less elaborate round entrance may be 
sufficient. Binnie(1938) was able to show that the use 
of a short rounded entrance was sufficient to suppress 
Orifice flows 

The use of an enlarged sloping entrance as adapted by 
the City of Edmonton, may be effective in suppressing 
orifice flow. Due to the lack of experimental 
observation, we may assume that if the lower nappe is 
able COvmhiil ithe.) pipeliativa, head..ftequal to the 
head-chamber height, orifice flow will not occur. Some 
experimental nappe profiles over ae sharp-crested 
circular weir are shown in Fig. 21. However, this inlet 
is not as efficient as the one described in (1) and it 
is not known whether cavitation will be a problem in 
the entrance. 

Binnie(1938), in his experiments, demonstrated that for 
a sharp-crested, re-entrant tube having a length more 
than 24 times its diameter, orifice control may not 
Govern, or that it may “contsol the flow only 
momentarily before it changes back to full-pipe flow. 
This 'length to diameter' control seems to be a result 


of free surface entrainment in the pipe with the air 
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Figure 21 - Typical nappe surfaces for circular weirs (From 
Wagner(1956) ) 


del 


Space being evacuated by the flow, hence prohibiting 
orifice flow formation. With zero sill height and an 
enlarged entrance, this required ratio could be much 


smaller. 


3.2 Experiments on Horizontal Circular Weir 

To verify the derivation of the weir flow equation (92) 
shown in section 3.1.2.2, an experimental investigation was 
Carried out. Details of the apparatus and results are as 


follows. 


3.2.1 Description of Apparatus 

A sketch of the apparatus used is shown in Fig. 22. The 
tank A, in which the outlet pipe D was clamped, wasS square 
in plan, three feet by three feet, and three feet four 
inches in height. The outlet pipe with four inch inside 
diameter was four feet four and a half inches long made of 
plexiglass. It was inserted vertically through the bottom of 
the tank and clamped approximately at the middle of the 
pipe. A circular G.I. sheet of twenty-three inches in 
diameter was laid concentrically on top of the vertical pipe 
D, with a hole four inches in diameter at the center. Above 
the circular plate there were four rectangular plates held 
vertical by two blocks of wood which partitioned the upper 
chamber into four equal regions. The four vertical sheets in 
the upper chamber stabilized the approach flow and 


eliminated any vortexes so that the discharge over the tube 
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Vertical discharge pipe D 


Inflow Pipe 


Circular plate(23"D) 


Partition wall 


Tank A 


(a)Plan view 


(b)Sectional view 


Diffuser 


from sump 


Figure 22 - A sketch of the circular weir apparatus 
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was radial. 

The inflow was pumped from a sump below the apparatus 
to be discharged into the lower chamber of the tank A. The 
inflow line was equipped with a bypass to control the _ flow 
rate. The end of the inlet was mounted with a diffuser which 
assisted in suppressing vortexes and turbulence. The 
discharge from the vertical pipe D returned to the same sump 
where water was being pumped. The inflow was measured using 
an elbow meter which was calibrated in situ. The calibration 
curve is shown in Fig. 23. The depth of flow at various 
locations above the circular plate was measured uSing a 


point gage. 


3.2.2 Experimental Observations and Data 

As the head rose above the crest of the pipe D, the 
flow at first adhered to the vertical pipe wall. It was 
clear and quiet with no air entrainment. A further increase 
in head caused the lower nappe to leave the pipe wall and 
close the center of the pipe. At the center of the pipe, a 
boil was formed above the upper nappe. The flow remained 
clear and quiet near the crown of the pipe which expanded 
and turned white and foamy as it approached the outlet. An 
increase in discharge also caused the boil to rise. The top 
of the boil must be kept below the end depth of the circular 
weir if the theory derived in this report is to remain 


valid. 
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Six sets of data were collected and tabulated in 
Table 2. The experimental results together with the 
theoretical weir flow equation were plotted in Fig. 24. If 
the boil rises above the critical depth of the approach 
channel, orifice flow will occur whose flow characteristics 
will be completely different from the weir flow discussed 


above. 


3.3 Design and Balance Sizing of the Inlet Sewer and _ the 
Drill Drop Pipe 
The-crvterva forgsizing drijdi- drop structures may be 
listed as follows: 

Ts Inlet size for the design discharge can be determined 
using the terminal-weir equation as discussed earlier. 

Zs The drop shaft should be designed in such a way that 
the design discharge will occur in the weir flow branch 
of the rating curve near the transition to full pipe 
a eo) 

Ore The inlet entrance should be designed to avoid orifice 
control. Also, the length to diameter ratio must 
satisfy conditions to limit cavitation. 

4, The head chamber depth should be at least equal to the 
head at design discharge, so that back-up will not 
occur at the inlet sewer. It may be desirable to keep 
the height larger than H as a safety factor against 
surcharge situation. 


oF The head chamber should be provided with an antivortex 
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Table 2 


Data and computation result for circular weir experiments 


(m?/s) (m?/s) 


Q(calc.) Q(measured) 
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device so that discharge through the drop pipe is not 
reduced by vortex formation. Two curtain walls across 
the diameter of the chamber will be sufficient to 


inhibit vortex formation. 


3.4 Air Entrainment 

Air entrainment of a circular wall jet in ae vertical 
drop pipe penetrating a free water surface was summarized by 
Rao(1975) which seems to support a relation of the form B = 
kF?, The experimental data is given in Fig. 25. In the 
relation, k was found to vary for different jet diameters. 
As explained in chapter 2, scale effect is inherent in most 
air entrainment model studies. 

The discussion of air entrainment for the shaft drop 
Structure is equally applicable to drill drop structure, 
with the exception of particulars due to geometric 


differences. 


3.5 Modifications and Recommendations for Further Studies 
Based on the above discussion, the following 

recommendations can be suggested: 

Ws Drop pipe entrance should be modified to avoid orifice 
control. Some approximate criteria to inhibit orifice 
control is given in section 3.1.4. Some experimental 
observations on the entrance geometry and the drop pipe 
Sizing is needed to confirm satisfactory operation. 


aie An anti-vortex wall should be provided to avoid vortex 
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Figure 25 - Air entrainment in a ringshaped wall jet (From 
Rao and KobuS( 1975) ) 


ie) 


Pas 
ere oie 
ave ee 


0 % | 
* ya a | ° 
~ ni : 
ie 
a | Le 
vy be 7 i 7; 
wae = os : 


80 


formation in the chamber. 

The entrance loss coefficient should be determined from 
experiments to facilitate designs. 

The appropriate chamber size to drop pipe size should 
be determined to avoid boundary effects. 

It was observed that the latter part of the weir flow 
curve changes into air-controlled flow. This type of 
flow iS geometrically specific and requires a model 


Study to accurately determine its rating curve. 
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4, Summary 

The study evaluates the two types of drop. structures 
used in the City of Edmonton. 

For the shaft drop structures, the inlet capacity is 
evaluated using terminal weir equations for circular 
channel. The entrance condition to the shaft is important in 
determining the types of flow in the drop shaft, or whether 
Cavitation is going to occur at the inlet, hence a smooth 
elbow is suggested in place of a 45° inclined floor. 

The flow of air-water mixture is analysed and suggested 
as one of the possible flow types in the drop shafts. 
Aeration of flow is a major concern in the analysis of 
enclosed drop structures, since it is the most important 
element controlling the capacity of the structure. Some 
aspects of air entrainment has been discussed. The analysis 
of two-phase phenomenon dynamically is extremely complicated 
and Gatrscalts even experimental approach have shown 
Significant scale errors. Therefore, it is suggested that 
further work be done in this area. 

The plunge pool is evaluated for effectiveness in the 
dissipation of kinetic energy before the flow entering the 
outlet. It is found that the plunge pool can be design aS a 
hydraulic jump basin or a jet diffusion basin. Suggestions 
for the geometry of the basin is also included. 

The effect of outlet submergence on the overall 
capacity of the structure, the problem with uplift pressure, 


overloading, balance sizing of various components are also 


81 


di ed ' ew, 
geub -jo 3 ows sag serait eva ‘ila. a7 
ra 72: ea 


“nogbe Ss +6 ody © ae 


, 4¢ - * © 
lines tra ovis qhab Siede:ef 
5 Bf 
a ' ~| 
ia ct 57i5 4 ; t Pd [23 BI a i 
: or ' 
P o- a = oe + oe ar 
i tiBa ta < OPT SS 20 DP etl Aw hb Ory, 
i 
| 
~ / i a e 7 = 
4 ~ a" ee ae ie | 19) ue 
— 
. ~ ‘ — or 
t 5 a) at 5 \ j iJ fe 4 
‘ on? ry 
" . — 
» 3 j tt) Sed24poge aeowedae 
a : ri 
= | 12 q 
E : ’ a as 7 a 
eixeroua bans ine (ih ta¢6e-Yis Jo wold gar ie 
Je a 4 
4 a H - > 
«27 T0n8 i ole = sidteeog a0 
- 93 i ¢ ? ‘“ ' 4 “x -~ i q >, 
, ‘ . 1 rh es | >” Bal a's a» 
i 
tral Voc {2 S. : aC. Reto oir ste qotD 
ome ~ : ; ot eee . ae , ae 
moe .% i? ee ‘T'o6den seid ber lloiwsne® 4 
- —g'j ; , 


aiayiene 2iT ,Gercioeih geen usa Ire: si9na tie 26 sApagae 
. ~~. = . ; ~~ a Ty 
o besa: £108GQ LeOROTT RS Bs Cees SiLMBAVG lonsnins nada ‘ikea 


nvoda evéed hac ie Ieatetstequs Oe ime PEA be me 
7 ¢ Fi i, nS ; iy a 7) 
gad? Sasteggue :.. 3 ‘e Sonat? paeyas om 2a “tewolt bagi 


4 Geis > isk ai S008 oe, a0 tedzay’ a 


Lvs 


efi i Z28¢4vii2stle 16) be nuda i ‘iia ial a. 


A ed 


ws 


ais en zane wollt ef4. -ssotad ter9ne stdwings ie. he 


ey . y= 


s sal npiess “ fas fooq snaydig eng ‘aul bys? Ree a 


: 7. nehseegtes, Miaaec ncie autite at al we ntwad ot 
7 ' "A 7 
on bat ton 73 oats. a neatid 0: , ~ iy 


- 7 
a) entedseaden ta o ~*~ 

| ) P . a 
rit 


- 
iv 


82 


discussed. 

For the drill drop structures, the circular weir flow 
equation is derived and verified with experiments. Weir 
flow is one of the flow types ‘that can exist in these 
structures. The other possible types of flow are the 
fuli-pipe: tit0w,, the cavitating flow, and the orifice flow. 
Orifice flow can be avoided by providing a smooth inlet 
entrance. Cavitating flow can cause excessive wear to the 
Pipe and should be avoided by proper design or by the 
provision of air vents. 

The various types of flow controls, their governing 
equations, and the transitions from one type of flow to the 
next are Summarized in a computer program. Discussion of air 
entrianment, balance sizing, and the control of orifice flow 
are also included. 

Based on the results of the Study, various 
modifications to the design of these structures have been 
Suggested. Recommendations for further studies, especially 
the air entrianment aspects and the energy dissipations 


aspects of these structures, have been provided. 
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Appendix I 


Dimensionless quantities for circular channels 
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Appendix II - Binnie's Description of Flow over a Re-entrant 


Tube 


Types of flow over a vertical re-entrant tube 


Binnie(1938) investigated the types of flow over a 


two-feet-long tube with one inch inside diameter. The tube 


was clamped to the base of a large tank at the middle of the 


Pipe. 


From his experiments, he identified the following 


possible flow types: 


1. 


Weir flow: As the head rises above the crown of the 
tube, water flows down the surface forming an annular 
jet with a glassy appearance. As discharge increases, 
the flow becomes more violent, of the noisy "gulping" 
type, with air entrainment at the center of the tube. 
Intermittent vortex may also form. The functional 
relationship is generally written as Q = f,(H?’?). 
Borda free flow: This is the same as_ orifice flow 
described in the main text of this report. This flow 
stage is quiescent with little air entrainment. The 
water surface in the tank is about level. As 
Binnie(1938) described, "The issuing stream assumed the 
appearance of a _ thin fixed glass rod projecting from 
the center of the pipe but not touching". 

Borda full flow: As the stage of the Borda free flow 
drops, the water surface above the tube suddenly opens 
up with the downward stream deflected to adhere to the 
tube surface. This iS Similar to the early stage of 


weir flow except for a small ring of air entrapped 
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around the inner surface of the top of the pipe. 

Full pipe flow: This occurs when the head above the 
tube becomes very high. The (tubeysflows full, quiet, 
with no air entrainment. This flow may exist at the 
Same stage as the orifice flow and which governs will 
depend largely on the geometry of the tube inlet. Some 
methods of preventing orifice flow are discussed in the 
main text. The functional relation of this flow is 


generally written as Q = f,(H'’?). 
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Appendix III - Computer Program for the Design of Drill Drop 
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Pipe and Sample Output 


THIS PROGAM UTILIZES THE HYDRAULIC CONCEPTS PRESENTED 
IN THE REPORT " A PRELIMINARY STUDY OF DROP STRUCTURES 
IN THE EDMONTON STORM WATER SYSTEM " TO PROVIDE A TOOL 
FOR THE ANALYSIS AND SIZING OF THE DRILL DROP STRUCTURES. 


MAIN PROGRAM 


COMMON DIAM, ENDPH(20), QDROP(20), AINIYE, GRAVTY, CHDIAM 
COMMON /SUBO2/ HWEIR(20), HFULL(20), HCAVI(20) 

COMMON /SUBO3/ALENGT,FPRIME, PRDOWN, ALOSSC(20),ALOSSB,PRCAV 
COMMON /SUBO4/ AMANN, IUNIT, QMAX, QMIN, ISTART 


DIAM = DIAMETER OF DROP PIPE 

ENDPH = END-DEPTHS 

Q@DROP = DISCHARGES 

AINIYE = INITIAL VALUE OF END-DEPTH 

GRAVTY = GRAVITATIONAL ACCELERATION 

HWEIR = HEAD FOR WEIR FLOW 

HEUE = HEAD FOR FUEL-PIPE FEOW 

HCAVI = HEAD FOR CAVITATING FLOW 

ALENGT = LENGTH OF DROP PIPE 

FPRIME = FRICTION FACTOR OF PIPE 

PRDOWN = PIEZOMETRIC PRESSURE HEAD AT OUTLET 

ALOSSC = ENTRANCE LOSS OF DROP 

ALOSSB = BEND LOSS OF DROP 

PRCAV = PRESSURE HEAD AT CAVITATING REGION 

AMANN = MANNINGS N 

QMAX = MAXIMUM DISCHARGE 

QMIN = MINIMUM DISCHARGE 

IUNIT = TYPE OF UNIT USED (BRITISH OR METRIC) 

ISTART = FLAG INTICATING TERMINATION OR GOING AHEAD 
CALL READ 


IF (ISTART .GE. 1) GO TO 99 
CALL DETRMN 
CALL PREPAR 
CALL PRINT 
CALL ENDPTH 
CALL HEAD 
CALL RESULT 
CALL PRINT2 
GO TO 1 
STOP 
END 
THIS SUBROUTINE WILL INTERPOLATE FOR THE CONTRACTION 


COEFFICIENT (ALOSSC). 


SUBROUTINE PREPAR 

COMMON DIAM, ENDPH(20) , QDROP(20), AINIYE, GRAVTY, CHDIAM 
COMMON. /SUBO3/ALENGT,FPRIME, PRDOWN, ALOSSC(20), ALOSSB, PRCAV 
DIMENSION VEL(14), 02D1(14), TABLE( 14,14) 

DATA MEU/O. 02.) 3s, 44 Sor eetere aOl 2a ton 20 1, 80n, 4067, 

HATA D204/7.0, 1. 4.9525 1, 4146, 18h 2-0, 202.2.0,9-0,4.0,9.0, 
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PIE1 = 3.14159265 
READ(3,501) ((TABLE(I,J),JU=1,14),1=1, 14) 
FORMAT (7F4.2) 


DO 550 K=1,20 
VELCTY = QDROP(K)/PIE1/((DIAM/2.)**2) 
IF (VELCTY.GE.40.) GO TO 561 
KK = 1 
IF ( VELCTY.LT.VEL(KK) ) GO TO 503 
KK = KK+14 
GO TO 502 
R3LATE = VEL(KK) 
R3ERLY = VEL(KK-1) 
JJ = 1 
RD2D1 CHDIAM/DIAM 
IF ( RD2D1.LT.D2D1(UJU)) GO TO 505 
JJ = JJ+1 
GO TO 504 
R4LATE = D2D1(uJu) 
R4ERLY = D2D1(JuJ-1) 
RANGE1 = TABLE(JUJ-1,KK) - TABLE(JUJU-1,KK-1) 
RANGE2 = TABLE(JJU,KK) - TABLE(JUJU,KK-1) 
CINTER = (VELCTY-R3ERLY)/(R3LATE-R3ERLY) 
FINTER = CINTER*RANGE1 + TABLE(JUU-1,KK-1) 
SINTER = CINTER*RANGE2 + TABLE(JUJU,KK-1) 
RANGES = SINTER-FINTER 
DINTER = (RD2D1-R4ERLY)/(R4LATE-R4ERLY) 
ALOSSC(K) = DINTER*RANGES + FINTER 
GO TO 550 
CNLOSS = 0.582 + 0.0418/(1.1-DIAM/CHDIAM) 


ALOSSC(K) = (1/CNLOSS-1)**2 
550 CONTINUE 


RETURN 
END 


THE FOLLOWING SUBROUTINE COMPUTES THE END-DEPTH FOR IN 
THE DRILL-DROP UPPER CHAMBER USING THE WEIR FLOW EQUATION 
DEVELOPED IN THE REPORT. 


SUBROUTINE 


ENDPTH 


COMMON DIAM, ENDPH(20), QDROP(20), AINIYE, GRAVTY, CHDIAM 
PIE = 3.14159265/12. 


ACONST 
BCONST 
YRATIO 
CCONST 
CONSTS 
CONST4 
CONST3 
YENDR 
xX 


SIN(PIE)**2 

2. *ACONST 

©On725 

2\3/ 4. *ACONST* (4. YRATIO IF VRATIOS* 2.0) 
16. * ( BCONST-CCONST ) 

8.*BCONST - 4.*CCONST 

BCONST 

AINIYE/(DIAM/2. ) 

VW2e 


THE FOLLOWING DO LOOP COMPUTES THE END-DEPTHS FOR 20 DIFFERENT 
DISCHARGES. NEWTON-RALPHSON TECHNIQUE IS USED TO SOLVE FOR 
YENDR. 
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DO 
KI 


100 
* 
* 
150 
200 co 
RETURN 
END 
(@: 
Cc 
GC DHE VEOEEG 
CG UPPER CHA 
C EQUATIONS 
C INPUT DIS 
Cc 
Cc 
SUBROU 
COMMON 
COMMON 
COMMON 
Bred 
DO 300 
C 
C WEIR FLOW 
C 
YEND 
QSTA 
HEAD 
HWE I 
Cc 
CG EUPESR LP E 
Cc 
VELC 
COEF 
VHEA 
HFUL 
Cc 
C CAVITATIN 
Cc 
HCAV 
300 CONTIN 
RETURN 
END 
Cc 
Cc 
C THIS SUBR 
C THE PARTI 
CYVSRESULTS.AF 
Cc 


200 I=1,20 
=) PG. 
DCONST =(QDROP(KI) **2/GRAVTY/(DIAM/2) **5/X**2) *YRATIO**2 
CONST1 = DCONST * 4. 
CONSTO = DCONST * 0.275 
FCN =(CONST3*YENDR**3+CONST4*YENDR**4+CONSTS*YENDR**5) 
-(CONST1*YENDR + CONSTO) 
DFCN =(3.*CONST3*YENDR**2 + 4.*CONST4*YENDR**3) 
+(5*CONST5*YENDR**4 - CONST1) 
YENDR2 = YENDR - FCN/DFCN 
DIFF = ABS( YENOR2-YENDR ) 
IF ( DIFF .LT. 0.0001 ) GO TO 150 
YENDR = YENDR2 
GO TO 100 
YENDR = YENDR2 
ENDPH(KI) = YENDR*(DIAM/2.) 
NTINUE 


WING SUBROUTINE COMPUTES THE HEADS IN THE DRILL-DROP 
MBER USING THE WEIR FLOW EQUATION, THE FULL-PIPE FLOW 
, AS WELL AS THE CAVITATING FLOW EQUATION FOR THE 
CHARGES. 


TINE HEAD 
DIAM, ENDPH(20) , QDROP(20), AINIYE, GRAVTY, CHDIAM 


/SUBO2/HWEIR(20), HFULL(20), HCAVI(20) 

/SUBO3/ALENGT,FPRIME, PRDOWN, ALOSSC(20), ALOSSB, PRCAV 
= 3.14159265 

I=1,20 


R = ENDPH(I)/(DIAM/2. ) 

R2 = QDROP(1)**2/GRAVTY/(DIAM/2) **5 
R = YENDR + QSTAR/8/PIE1**2/YENDR**2 
R(1) = HEADR*DIAM/2 


FLOW 


TY = QDROP(I)/PIE1/(DIAM/2.)**2 

F = 1 + FPRIME*ALENGT/DIAM + ALOSSC(I) + ALOSSB 
D = VELCTY**2/2./GRAVTY 

L(I) = COEFF*VHEAD + PRDOWN - ALENGT 


G FLOW 


I1(I)} = PRCAV + (ALOSSC(I) + 1)*VHEAD 
UE 


OUTINE DETERMINE WHICH FLOW TYPE GOVERNS FOR 
CULAR DROP PIPE SIZE AND DISCHARGE. 
OR THE RUN WILL ALSO BE PRINTED. 
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SUBROUTINE RESULT 

COMMON DIAM, ENDPH(20) , QDROP(20), AINIYE, GRAVTY, CHDIAM 
COMMON /SUBO2/HWEIR(20), HFULL(20), HCAVI(20) 

COMMON /SUBO3/ALENGT,FPRIME, PRDOWN, ALOSSC(20), ALOSSB, PRCAV 


DETERMINE THE GOVERNING FLOW TYPE 


405 
415 


416 


420 


417 


410 


418 
400 


DO 400 I=1,20 
IF (HWEIR(I).GE.HFULL(I)) GO TO 405 
IF (HFULL(I).GE.HCAVI(I)) GO TO 410 
GO TO 415 
IF (HWEIR(I).GE.HCAVI(I)) GO TO 420 
HACTUL = HCAVI(I) 
WRITE (6,416) QDROP(I), HACTUL, ALOSSC(I) 
SORMATEN “O° «10X,.FS.2, 13K bole Tike CAVITATING’, 10X,F 12.3) 
GO TO 400 
HACTUL = HWEIR(I) 
WRITE (6,417) QDROP(I), HACTUL 
FORMAT (207 3 10X,FS8n20 (aX bes texec WETR FLOW”, 20X,°--") 
GO TO 400 
HACTUL = HFULL(T) 
WRITE (6,418) QDROP(I), HACTUL, ALOSSC(T) 
EFORMAIN CO" . 10X.FS.2.13X%.F 7.2. 12%, FULL. FLOW”, 10X,F 12-3) 
CONTINUE 
RETURN 
END 


THIS SUBROUTINE READS IN THE NECESSARY DATA INTERACTIVELY 


101 


103 


104 


105 


106 


107 


108 


102 


109 


110 


SUBROUTINE READ 

COMMON DIAM, ENDPH(20), QDROP(20), AINIYE, GRAVTY, CHDIAM 

COMMON /SUBO3/ALENGT,FPRIME, PRDOWN, ALOSSC(20), ALOSSB, PRCAV 

COMMON /SUBO4/AMANN, IUNIT, QMAX, QMIN, ISTART 

WRITES C7104) 

FORMAT (’THIS IS A PROGRAM DESIGNED TO ANALYZE DRILL DROP’ ,/ 
*’SIZE AND CAPACITY. IT IS DESIGNED TO HAVE ALL THE DATA’,/ 
*’ INPUTTED INTERACTIVELY. ’) 

WRITE (7,103) 

FORMAT (‘ENTER DROP PIPE SIZE, LENGTH AND ROUGHNESS’ ,/ 
*/IN THE FORM 2*XXX.XX, X.XXXXX’,/ 

*/OR HIT BREAK TO TERMINATE THE RUN’) 

READ (7,104) DIAM, ALENGT, AMANN 

FORMAT (2F6.2, F7.5) 

WRITE (7,105) 

FORMAT (‘ENTER HEAD CHAMBER DIAMETER, XXX.XX’) 

READ (7,106) CHDIAM 

FORMAT (F6.2) 

WRITE (7,107) 

FORMAT (’ENTER BEND LOSS COEFFICIENT OR ZERO, XX.XXX’) 

READ (7,108) ALOSSB 

FORMAT (F6.3) 

WRITE (7,102) 

FORMAT (‘ENTER OUTLET PRESSURE HEAD, XXX.XX‘) 

READ (7,109) PRDOWN 

FORMAT (F6.2) 

WRITE (7,110) 

FORMAT (’ENTER MAXIMUM AND MINIMUM DISCHARGE, DEXXXX XX” ) 
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READ (7,111) QMAX, QMIN 

FORMAT (2F7.2) 

WRITE (7,112) 

FORMAT (’ENTER INITIAL VALUE OF END DEPTH XXX.XX’) 
READ (7,113) AINIYE 

FORMAT (F6.2) 

WRITE (7,114) 

FORMAT (‘’WHAT UNIT ARE YOU USING?’/ ‘ENTER O FOR BRITISH, 
*4{ FOR METRIC’) 

READ (7,115) IUNIT 

FORMAT (11) 

WRITE (7,116) 

FORMAT (‘DO YOU WANT TO GO THROUGH WITH IT?’/ 
*/ENTER O FOR GO AHEAD, 1 TO NEGATE’) 

READ (7,117) ISTART 

FORMAT (11) 

RETURN 

END 


THIS SUBROUTINE DETERMINE THE TYPE OF UNIT TO BESUSEDE 

THE ARRAY OF DICHARGES WITHIN THE RANGE ENTERED, THE CAVITATING 
PRESSURE, ALSO IT WILL CONVERT MANNING’S N INTO DARCY-WEISBACH 
FRICTION FACTOR. 


700 


710 


720 


SUBROUTINE DETRMN 

COMMON DIAM, ENOPH(20), QDROP(20), AINIYE, GRAVTY, CHDIAM 
COMMON /SUBO3/ALENGT,FPRIME, PRDOWN, ALOSSC(20), ALOSSB, PRCAV 
COMMON /SUBO4/AMANN, IUNIT, QMAX, QMIN, ISTART 

QINTVL = (QMAX-QMIN)/20. 

QDROP(1) = QMIN 

DO 700 I=2,20 

QDROP(I) = QDROP(I-1) + QINTVL 

CONTINUE 

TF CYORET {FO . 1) GO TOTO 

GRAViIiVea=n so 22 


PRCAV = -31. 

GO TOFT720 

GRAVTY = 9.81 

PRCAV = -9.45 

CONT INUE 

FPRIME = (185.*AMANN**2)/(DIAM**(1./3.)) 
RETURN 

END 


THIS SUBROUTINE WILL ECHO-CHECK INPUT DATA AND PRINT HEADINGS 
FOR THE RESULTS. 


808 


801 


SUBROUTINE PRINT 

COMMON DIAM, ENOPH(20), QDROP(20), AINIYE, GRAVTY, CHDIAM 
COMMON /SUBO3/ALENGT,FPRIME, PRDOWN, ALOSSC(20), ALOSSB, PRCAV 
COMMON /SUBO4/AMANN, IUNIT, QMAX, QMIN, ISTART 

WRITE (6,808) 

BURMATEG IS. f ffi) 

WRITE (6,801) DIAM, PRCAV 

FORMAT (10%, 9X, “DIAMETER OFSPIPE, =", He 26 MOS 

*’CAVITATING PRESSURE =’, F8.2) 
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300 WRITE (6,802) ALENGT, ALOSSB 

301 802 FORMAT (‘0’, 9X, ‘LENGTH OF DROP =’, F8.2,10X, 
302 *’BEND LOSS COEFF =z’, F8.2) 

303 WRITE (6,803) AMANN, PRDOWN 

304 803 FORMAT (’0’, 9X, ‘MANNINGS N =’, F8.4,10X, 
305 */OQUTLET PEIZ. HEAD =’, F8.4) 

306 IF (IUNIT .EQ. 1) GO TO 850 

307 WRITE (6,804) CHDIAM 

308 804 FORMAT (’0’, 9X, ‘HEAD CHAMBER SIZE =’, F8.2,10X, 
309 *’UNIT IS BRITISH’) 

310 GO TO 860 

311 850 WRITE (6,805) CHDIAM 

312 805 FORMAT (‘0O’, 9X, ‘HEAD CHAMBER SIZE =’, F8.2,10X, 
313 *’UNIT IS METRIC’) 

314 860 WRITE (6,806) 

315 806 FORMAT (’0’,8X,78 (’*’)) 

316 WRITE (6,807) 

317 807 FORMAT (’0’, 10X, ‘DISCHARGE’, 15X, ‘HEAD’, 10X, 
318 *’ TYPE OF FLOW’, 7X,’ENTRANCE LOSS COEFF’) 

319 WRITE (6,806) 

320 RETURN 

321 END 

322 c 

323 c 

324 C THIS SUBROUTINE PRINTS THE HEADS FOR THE INPUT DISCHARGES 
325 C FOR WEIR FLOWS, FULL-PIPE FLOW , AND CAVITATING FLOW. 
326 Cc 

327 € 

328 SUBROUTINE PRINT2 

329 COMMON DIAM, ENDPH(20) , QDROP(20), AINIYE, GRAVTY, CHDIAM 
330 COMMON /SUBO2/HWEIR(20), HFULL(20), HCAVI(20) 

331 WRITE (6,912) 

332 WRITE (6,911) 

333 WRITE (6,910) :; 

334 910 FORMAT (’0’,8X, ‘DISCHARGE’,10X,’WEIR HEAD’,7X, 
335 *’FULL-FLOW HEAD’ ,7X, ’CAVITATING HEAD’ ) 

336 WRITE (6,911) 

337 911 FORMAT (’0’,8X.71('*’)) 

338 912) FORMAT €74"../?=2 7°77) 

339 60 °(S20 1 2177 20 

340 WRITE (6,913) QDROP(I), HWEIR(I), HFULL(I), HCAVI(I) 
341 913. FORMAT (‘0’,8X,F9.2, 10X,F9.2,7X,F14.2,7X,F 15.2) 
342 920 CONTINUE 

343 RETURN 

344 END 


End of file 
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DIAMETER OF PIPE = 1.00 CAVITATING PRESSURE = -31.00 
LENGTH OF DROP = 100.00 BEND LOSS COEFF = 0.0 
MANNINGS N =O O30) OUTLET PEIZ. HEAD = 0.0 
HEAD CHAMBER SIZE = 4.00 UNIT IS BRITISH 
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100.00 1024.72 FULL FLOW 0.341 
95.00 915.06 FULL FLOW 0.341 
90.00 3h ql gO FULL FLOW Owsan 
85.00 Val2eion FULL FLOW OS 4:4 
80.00 CSesZ FULL FLOW 0.341 
1S) C18; 532.65 FULL FLOW 0.344 
70.00 ASiim tat FULL FLOW 0.341 
65.00 SYS). VS FULL FLOW 0.341 
60.00 304.90 FULL FLOW 0.341 
55.00 240.23 FULL FLOW 0.341 
50.00 esta} ste} FULL FLOW 0.341 
45.00 VOT a Ih FULL FLOW OE 
40.00 79.95 FULL FLOW 0.341 
35.00 Sie FULL FLOW 0.341 
30.00 Wah). FULL FLOW O}. SENS: 
ZOOO 0.89 WEIR FLOW me 
20.00 0.78 WEIR FLOW = 
VSIOe) 0.66 WEIR FLOW = 
10.00 OFS2 WEIR FLOW ies 
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